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Needs [ "Laurent Functions * H ] 



Needs [ w Laurent Not at ionTes t * " ] 

Needs: :nocont : Context LaurentNotationTest • was not created when Needs was evaluated. 
We develop the modulating function that is used in GSM first to use as an example to demonstrate Laurent's idea. 



3 



812500 
BT := 0.3 



Modulationlndex : = — 
2 



See LaurentTheory.nb to see the derivation of the C functions 

When running the notebook we fix the value of L here and the rest of the graphs and functions use this value 



Office vU; tftte.-U^X^ d if^.V A w OOCiA-S \ \ l*L-«.,Z\k-L., Ac,)c 







J> A .- /J^jr 
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£ : = 6; 

PhaseFunction[t_] = N[0i,,t]; m . „ , 

phasepoints = Table [{t, Pha»eFunction[tT]} //*. {t, 0, I., 1/40)]; 
l.istPlot[pnasepoints, Plot Joined -> True] 

"Sfgra^ralLa to converge to prescribed accuracy after 7 recursive bisections 
in LaurentFunctions' Private • tl near LaurentFunctions ' Private tl = -1.91796x10 . 

N £nte g ra?e : ?allea to converge to prescribed accuracy after ' -cursive bi^tion. 
in LurentFunctions' Private >tl near LaurentFunctions Private tl - -5.52869x10 . 

integrate 1 railed to converge to prescribed accuracy after 7 recursive bi.-cti««. 
in LaurentFunctions' Private <tl near LaurentFunctions -Private < tl = -1.96696x10 . 

further ! output of NIntegrate: :ncvb will be suppressed during this calculation. 

integrate: =slwcon.: ^^ft^^tSS^ i^uf ft^ent^^glrel^sfon^' ^ 
% S^ 1 ?5SK5irii"c2imrw t S? usinf the' option M ethod->Oscillatory in integrate. 

value 




If your integrand 



NIntegrate. 



G |CrSer ! output of NIntegrate: : slwcon will be suppressed during this calculation. 




1 2 3 4 5 6 

- Graphics - 

Bitseq = Table [1, {i, 1, 20}] 

{1, 1, 1. 1. 1, 1- 1- 1- I- 1. !< X ' X ' X ' 1 - X ' X ' X ' 1} 

RandomBitSea; = Table [Random [Integer, {0, 1}], {i, 1, -40}] // Map[# < -2) + 16, #]* 



{ 



ii-i -l. -l, i, i. i- !< 1 - 1 - - 1 ' 1 - V -,V 11 _1, _1, _1 ' 
-i, i. i! -I. -i! -1. -i. 1. 1. 1 - - 1 ' - 1 ' - 1 ' 1} 



Testing The Modulator ... 

Calculating the bandwidth of the signal 

xll = Table[LaurentC[L][0][tT], {t, 0, X. + 1 - 1/200, 1/200}]; 

((Drop[xll, -1] -Drop[xll, 1]) 200 ) *2 / 200 // Apply [Plus. #]* 

1.29684 

Here the pulse 6T is just stretched to 8T to see the performance 



LaurentModulatorPulses.nb 



£xample of stretching a C 0 pulse Of 6 to 8 
In this section we consider filtering the pulse 

fil =,BesselFilter[8] [LowPass3dB [1 / t 2 Pi 0 . 3] ] [FrequencyResponse] W , 

Plot [20 Log [10, Evaluate [ (f il /. s 2 Pi X £ / T W / ah« n /<= n 
GridLines -> Automatic/ kxesLabel -> { - l/T^ei^ ] 3 ' ' ' °' 2> ' 

dB 

1/T 



-40 




Graphics 



Res[t_] = InverseLaplaceTransfonnffil, s, t] 



2027025 ( (9.36812 + 0. I) e~ 897172 * fc Cos[139301 tl 

(13.7944 .0. I) E-""«-'cos[420044. t] + (4.73111.0 I) E -"i35B. t Cos 
(0.304883 + 0. I) E-«"».t coe[l. 02016x10* tl + Cos[708768. t 

(29.1264 + 0. I) £-897172. t gj[ n [139301 t] - 

Plot [Res [ t T] , {t, 0, 6}, PlotRancje -> All] 




- Graphics - 



filTaps = Table [Res [tl], {t, 0, 6, 1/20}] //Chop; 
We estimate the group delay to be 1.8T 



Length [ f i lTaps J 
121 
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L := 8; 

Pb.aseFunction[t_] . = M[fa.tJ» 1/4011 . 
phasepointa = Table [{t, PhaseFunction[t T] } // N. {t, 0, I. , 1/40}]. 
I.istPlot[phasepoints, Plot Joined -> True] 

N ^t!g"te : ?aUek to converge to prescribed accuracy after 7 recursive W«ctl««. 
in LurentFunctions • Private -tl near LaurentFunctxons -Private tl = -3.68196x10 . 

^Igrate^allea to converge to prescribed accuracy after 7 r-ajr-iv- 

in LurentFunctions- Private -tl near LaurentFunctions Private tl - -7.4914x10 

integrate 1 falLd to converge to prescribed accuracy after 7 "c-sive bi-.eti««. 
in LurentFunctions -Private -tl near LaurentFunctions Private tl - -5.43232x10 . 

G Further "utput of NIntegrate: :ncvb will be suppressed during this calculation, 
integrate: =slwcon : Numerical integration converging too slowly; ^J^^^' ^ 
%\ ^Jt^^^<^^^ u fing r tne'o P r tion^^d->Oscillato r y in HXntegrate. 

integrate: :sl„con. : ^°££^% t <£^™ L%°uf f ^ent^oSgPrelXston^' ^ 
If your^egrand^s^scU^atory^ uslng'tnt option M ethod->Oscillatory in NIntegrate. 

NIntegrate: :slwcon.= ^^ffi^S^ i-uff i^^nt^&lrecTsf c^' ^ 
If ^^SSnTii^iSStS^ using^nf' option „ethod->Oscillatory in NIntegrate. 

G |urther : ou?pu : t of NIntegrate: : slwcon will be suppressed during this calculation. 
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1 
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2 4 6 B 

- Graphics - 

pulse = Table[LaurentC[8] [0] [t T] , {t, 0, 9, 1/20}]; 
pulse2 =Table[I.aurentC[8] [1] [tl], {t, 0, 7, 1/20}]; 
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FiltereaPulse - T/20 MapConvolve fpulse. fil Taps] //Take[# . {3 4, l.en gth[#] )J& 

2.10426xl0" 9 , 4 38922 x 10-* R Qqco^ , n-9 '*'* iyb8xl 0 , 9 . 86455 x 10 10 , 

i ,: M , , j;.S-: :.:»"S."^: ^I^v j-isj;. 1 ;^ «■«>»«■.»-. 

0.0029439, 0.00390362; 0 OOSImII C > Oolkwl**^^^ 0162196 ' 0.00219697, 
0.0172954, 0.021476, 0 '. 6264451 0 i 032~ 30of 0 I 57 n 06 „ 8 ;,S; 0109288 ' 0.0138093, 

0.0665603, 0.0783079, 0.0914948 6 lOfilfls" S'?i?}$° 7 n' °;2 470666 ' °-Q561754, 
0.20329, 0.22737, 0.25287 C I 279682 n Vn«Si "n*, 3 ,', 0 - 140261 ' 0.159691, 0.180711 
0.427934, 0.459013. 0.490013 0 520668 0 S^4n°p 3 ^ 6 f lli^ 3 6653 9 ' 0.397029. ' 
0.65922, 0.682107, 0 702805 O ^fnoB n -it^I 08 ' 0.579856, 0.607837, 0.634379 
0.770616, 0.771369 6 768969 o ?fiS2?4 VlfP' °- 7 49713. 0.759727, 0.76672 
0.692108 0.670178,' 0 S ' o"e!o 44 n^of a!? 6 «° " 743274 ' °-728852, 0 711736 
0.473326, 0.442277 0 Allot' S'SiSlls O «n«i Vf 429 ' °-534607, 0.504193.' 
0.239519, 0.214852 0.191636 0 llSSfl' O^fofo?' S' 32 f 221 ' °-29283, 0.26555 
0.0848118, 0.0723143 0 0612O1S o nJnoe,"? 7 ?!' 0.13122, 0.114217, 0.0987613 

§:§ol„^r„1t!™ Ili^^MsSAtp^"" 5 ' SiSSSSi^.-SfSSSfe,. 

! >«! »-««l»l"10 '. 1.4722x10-', 8.02333x10-". 4 27388x10-" 

=: o» ^lllUZ"' VSJL"^ —"x"-". i.iSii".-". 

,/. 52915x10 . 5. 81S84xlO 2S , 1. 75371xl0- 25 , 1.2 2 3 9 2X10" 26 01 
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FilteredPulse2 = T/20 MapConvolve[pulse2, filTapB] //Take[#, {50 ? Length!*] )]& 

/l 27972 xlO" 7 2 18306 x 10" 7 , 3 . 66972 x 10" 7 , 6.0806x10"', 9 . 93427 x 10" 7 , 1. 60076 x 10- 
{1.27972x10 . 2.1BJU6 x , « q 4505 9xl0-, 0.0000142688, 0.0000212804, 

0.00618453, 0.00739101, J^SSIm n n9AMl q 0272921 0 0297419, 0.0321652, 

S:!BSS: S:S : g |:S j.VS S j .««« . -.g".^. 
S:SSSI : S:StlfI||: : : » : S VKiiSJ- S-8ii!i;?o?o!Si!S!- . 

0.0306913, 0.0283283, 0.0259611, 0-0236189, 0 . 021328^ O uxsii ^ 0 . 00579518, 

0.0149823, 0.0130987, 0-0113504, Q ° -0097438, 0 • °° 8 282|| i " - 7 0.00105885, 

0.00476282, 0 00386365 . 0.0030899 0 ^""JJ^/^S'JgSi^figB, 0.0000870728, 
0.000765388, 0 . 000536896 0 . 000363597, 0 ' • £"£^23, 0 . 0000338422 , 0.0000428466, 
0.0000507313, 0.0000318118, 0-0000253258, 0-0000271231 o Q42 0 . 000 0692325, 

8-SSS5iI5ii5: S:SSSS Hit : SiSSSSixxSii: S:S?SS3SliS 5 : 0 0000330554 o o?oo236i 2a . 

0 : SSSSSJ"! 5 . 95924 x 10-; -1 - 64542 x 10^ - - 8 -^^7^ 

n l£lt«10-<' 7 16306x10-, 6.89804x10-, 6.49223x10-, 5.98586x10-, 

1 ttlllAo-'' I"5"mx10-« 4. 20863x10-, 3.6208x10-', 3.06728x10- 2.55958x10 , 
1^71xl2-' l'. 058 1. 36297X10-, 1.07376x10-. 8-34162x10- , 6.3908x10 7 , 
HVb\AI< 3.59854x10" 7 , 2.64486x10- 1 . 91715 x 10" , 1.370 44>« 0 • 

q «984 X 10- e 6 71324x10-, 4.59913x10-, 3.10539x10 , 2.06612x10 , l.Jb4ZXJ.u 
ll£ls*lo->' 5 55593x10-, 3.47564x10-, 2.13956x10- 1.29613x10-, 

i-SKS* fcSES=: £323^: 

IloioixlS- 17 ' 21B469K10-", 1.39012x10"", 1.92164x10" , 6.31946x10 » 
-6 "328 "o-». I'-^oLlO-, 1.42063x10-, 7 . 33335 x 10"» , 1.13717x10"", 

3.90929xlO" 19 , 1.30285xlCT 19 , 0} 

Int erpol at ion 

InterpolatingFunction[{{0. 0.0000492923}}, <>] 

""•SCR! Eiii < SS££i£^i A122i£24 / / ™ / / 

Interpolation 

Inter P olatingFunction[{{0, 0.0000389538}}, <>] 
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Plot [FiltPulse [ 8] [0] [t T] i {t, 0, 9}, PlotRange -> All] 
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2 4 6 8 

- Graphics - 

Plot [LaurentC [ 8 ] [0] [t T] f {t, 0, 9}, PlotRange -> All] 

0.8 
0.6 
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0.2 




- Graphics - 
Show[%, %%] 
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Graphics 



Plot [Fi It Pulse [8] [1] [t T] / {t, 0, 7}, PlotRange -> All] 



0.04 
0.03 
0.02 
0. 01 




1 2 3 4 5 6 7 
- Graphics - 

Plot [LaurentC [8] [1] [t T] / {t, 0, 7}, PlotRange -> All]. 

0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 




- Graphics - 

Bandwidth [TestPulse_, Supp_] : = 
Module [{xll), xll = Table [TestPulse [t T] , {t, 0, Supp - 1 / 200, 1/200}]; 
((Dropfxll, -1] -Drop[xll, 1]) 200 )*2 / 200 // Apply [Plus, #]&] 

Bandwidth [FiltPulse [8] [0] , 9] 

0.711659 

Bandwidth [LaurentC [8] [0], 9] 
1.29684 

Bandwidth [FiltPulse [8] [1] , 7] 

0.00338446 

torn = Modulator [L] [Table [1, {i, 1, 100}], Samplinglnterval -> T/10, 
NumberOf Curves -> 2, ModulatingPulse-* FiltPulse]; 

torn = Modulator [L] [RandomBit Seq, Samplinglnterval -» T/10, NumberOf Curves -> 2, 
ModulatingPulse -* FiltPulse] ; 
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ListPlot[Im[tom] , P lot Joined -> True] 
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- Graphics - 

I*istPlot[Re[toin] , Plot Joined -> True] 




D00 



- Graphics - 

riistPlot[Abs[tom] , PlotJoined -> True] 



0 00 




- Graphics - 



LaurentModulatorPulses.nb 



_10+) 



{Re [torn) , Im[tom] ) / / Transpose 



// ListPlot [#, PlotJoined -> True, AspectRatio -> 1]& 




- Graphics - 

torn = Modulator [L] [Randoms it So*, Samplinglnterval -> T/10, Number Of Curves -> 1, 
ModulatingPulse FiltPulse] ; 

ListPlot[Ixn[tom] , Plot Joined -> True] 
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- Graphics - 

ListPlot[Re[tom] , Plot Joined -> True] 



0.75 
0.5 
0.25 



-0.25 
-0.5 
-0.75 



)0 



600 



WWW 



00 




000 



- Graphics - 
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ListPlot [Abs [torn] , Plot Joined -> True] 




- Graphics - 

{Re [torn], Im[tom]} // Transpose // ListPlot [#, PlotJoined -> True, AspectRatio 




- Graphics - 

xll = Table[TestPulse[6] [0] [t T] , {t, 0, 7 - 1/200, 1/200)]; 

((Drop[xll, -1] - Drop [xll, 1]) 2 00 )*2 / 200 // Apply [Plus, #]& 

0.92632 

TestPulse[6] [0] [t_] = LaurentC[4] [0] [t 5/7 ] 

Tcnsr ♦■¥■]] ?skss + *f-\\ 

Sin[*[«. -LI]] 



8X2500 
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L := 5; 

PhaseFunction[t_] = 4>i*,t. * . 

phasepoints = Table [{t, PhaseFunction[t T] } // N, {t, 0, L , 1/40}]; 
ListPlot [phasepoints, Plot Joined -> True] 

NIntegrate: mcvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 1 . 48053 x 10" 6 . 

NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = -1.12255x10" . 

NIntegrate: mcvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 1 . 08808 x 10' 6 . 

further* output of NIntegrate: :ncvb will be suppressed during this calculation. 

NIntegrate- slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand, or insufficient WorkmgPrecision. 
If your integrand is oscillatory try using the option Method- >Oscillatory .an NIntegrate. 




Graphics 



Ii : = 6 

TestPulse2[6] [0] [t_] = LaurentC [4 ] [0] [t 6/7] 



1. Sin[#[4, 
Sin[*[4. 

torn = Modulator [L] [RandomBitSeq, Samplinglnterval T/40, NumberOf Curves -> 1, 
ModulatingPulse FiltPulse] ; 

Save["ModulatorData.m", {T, L, RandomBitSeq, FiltPulse, torn}] 
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I*i s t Plot [Im[ torn] , PlotJoined -> True] 
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- Graphics - 

ListPlot[Re[tom] , Plot Joined -> Tame] 
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- Graphics - 

ListPlot [Abs [torn] , PlotJoined -> True] 




00 



- Graphics - 
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{Re [torn], Im[tom]} // Transpose // Lis tPlot [#, PlotCToined -> True, AspectRatio -> 1]& 




- Graphics - 



xll = Table[TestPulse2[6] [0] [t T] , {t, 0, 7 - 1/200, 1/200}]; 
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Needs [ "Laurent Functions * ■ ] 

RuleDelayed::rhs : Pattern t_ appears on the right-hand side of rule 

PhaseAngle[L_j [t_] » (PhaseAngle [L] [t__] = Module [{xl, x2, x3 , x4 x5, x6} ( «1» ] ) . 

Needs [ " LaurentNot at ionTest * H ] 

Needs: rnocont : Context LaurentNotationTest • was not created when Needs was evaluated. 

Information on the functions used can be obtained using help. 
Names ["Laurent Functions 

{AKN, AlphaKI, ANKInitialStateSetUp, BT, FiltPulse, h, hFiltered, InitialState J 
SuS LaUren o!; K ' LaurentS, M, ModulatingPulse, Modulationlndex Adulator 
S^M^ r r S ', Ph c SeA 2 glS ^ Receiver , ReceiverProper, S, Samplinglnterval ' 
StartmgQuadrant, SyncSample, T, C, if/} 



T := 

812500 
BT := 0.3 



Modulationlndex : = i 
2 

<< ModulatorData.m; 



RandomBitSeq 

[1, 1, 
-1, 1, 



{1 ; h - 1 : -i- i- -i. i. -i. i. -i, i, -i, -i, -i, -i, i, i, i, -i, -i, i, lt lf 



7 -, '. A ' : A « -i« i» i» -i» -i, i, -i ( i, -i, if i/i,i. -l. -l -i i -i 

i, i, i, -i, i, -i. -i. i, x, i, .i, _x, 1( !, lt 1 - t 1 - i ^ 
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M-l 



With L = 8 and M = 2 we can utilise only the main values. It so happens that when K = 1, the dominant value occurs at 
2.5T in the filtered pulse. The three dominant values when K =0 occurs at 3.5T, 4.5T and 5.5T. The value of the pulse 
at 6.5T is smaller than the value of the pulse with K = 1 at 2.5 T. Thus the expression becomes 



§wt + at = J* 0 '"" 1 Pulse[0][3T + <5T] + J**-"-* Pulse[0][4T + 6T] + 
j*o.«-5 p u lso[0] [5T + ST] + J**-"-* Pulse[l] [2T + <5T] 

§ mt + at= J* 0 *"-' Pulse[0][3T + <5T] + J* 0 -"-* Pulse[0][4T + 6T] + 
j*o,n-s pulse[0] [5 T + <5T] + j^-"- a ••*->) Pulse[l] [2 T + <5T] 

Given that a N . 3 has been already decided, we can precalculate the possibilities and store them. 

g T= j*o,»-s { j{a H -4*o«i-i) Pulse[0][3T + <5T] + J 0 * -4 Pulse[0][4T + <5T] + 
" T *Pulse[0] [5 T + <5T] + J"*-** •"-» •«*-»••"-> > Pulse [1] [2 T + <5T] ) 

But J 2 = 1 and so we get 

§ HT + AT = jXo, H -5 ( * ♦««-»> Pulse[0] [3T + <5T] + J**"- 4 Pulse[0] [4 T + 6T] + 

HT + Pulse[0] [5 T + 6T] + tf«-«-4» > Pulse [1] [2 T + 6T] ) 

IiOokUpTable[Pulse_, <5T_] : = 
Module [{xl, x2, x3, x4}, 

x \Zi { . l -i?'-i}: { -.\ } ; f;i',; a. ^ i - «>» 

x2 [bltseq ] s = 

, jtoits^ausu ♦bttwCtaj] Pu i so [03t3T + <5T] + o»»*-*~«t I»3 1 Pulse[0][4T + 5T] + 
Pulse[0][5T + <5T] + j*it««iii]l [ij) pulse[l][2T + 6T] ) ; 

{xl, Map[x2[#]&, xl]} // Transpose] 

tab = LookUpTable[FiltPulse[8] , 0.5 T] 



fff-1 -1 -1) -0.00122183 - 0.770616 1}, {{-1, -1. 1}, 0.702339 + 0.770616 1}, 
ML!' 1 -1} 0.614125 - 0.770616 1}, {{-1, 1, 1}, 0 . 0869922 + 0 . 770616 I , 
{ 1, -l! -1 0 0869922 - 0.770616 1), {{1, -1, 1}, 0 . 614125 + 0 .770616 I 
{{l! 1,-1); 0.702339 - 0.770616 1), {{1, 1, 1}, -0 . 00122183 +. 0.770616 I) ) 

Sort [tab, #![[!/ 3]] >#2[[1, 3]]&] 




{{ 

FiltPulse[8] [0] [4 T + 0.5 T] 

0.770616 

Now we build a receiver! ! First generate the modulated sequence 
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Options [Modulator] 

{startingQuadrant -» 0, 
InitialState-» {1, 1, 1, 1, 1, l, 1, 1, i, i, i, 1# lt lt lf 1# lt 1# lf 1}/ 

SamplingInterval-> 2600^000 ' Number ° f Curves 4 , ModulatingPulse -* Laurentc) 



2600000 
1 

32 



00 / 



tom2 = Modulator [L] [RandomBitSeq, NuinberOf Curves -+ 2, ModulatingPulse -+ FiltPulse] ; 
ListPlot [{Re[tora2] , Ixn[tom2J) //Transpose, Plottfoined -> True] 

- Graphics - 
RandoxnBitSeq 

{1 i x i " x i x ; x ; ^ - 1 ' ^ - 1 ' - 1 - - 1 ' 1 > 1 > x - - 1 * 1, i, 1, 

; \ 1 ' r 1 ', 1 ' _1 ' x - 1 ' ^ -i- i. i. -i, -i, -i, i, -i, 

i, i, i, i, -i t -i, i, i, i, -i, -i, i, lt lt 1# lt lf lt lt lt _ lf _ lf 1# _ x ^ _ 1} 

Receiver[L] [tom2, StartingQuadrant -> 4] 

{1,1, 11 11-1-1-1 11, -1,1-1 l, -1, 1, -1, -l, -i, -x, i # lf 1# _ lt _ lt 

1, 1,1,-11-1 1 -1 1, -1, l, l, _i, -i, 1# _ 1# x _ x ! ! x i # . _! 

-1,-1, 1, -1, 1, 1, 1, -1, 1, -1, -l, l, !, i, _ a , 1# 1# l/lf 1# 1# lt ' lt lf _ 1} 

We now calculate the BER given the and L = 8. We write all the possible terms 

j*o,»-5 ( j(a B . 4 + ♦ Pul8e[0] [ 6T] + 

j^-^a.- , + ,*a..i) pui Be[0 ] [ T + <5TJ + j( a »-«* a w-i + pulse [0] [2 T + <5T] + 

j(«-.4*an-») Pulse[0][3T + <5T] + J«h-4 Pulse[0][4T + <5T] + Pulse[0][5T + <5T] + 
j(-a K - 5 ) pulse[0][6T + <5T] + j(-*»-s - ««-« ) Pulse[0][7T + <5T] + 

j(-a w .5-o H -« -a M . 7 ) puigejQ] [ 8T + 6t j + ^f( a n - 4 ♦ *h - j ♦ °k - j * on - i ♦ *s * *u - i ) Pul S 6 [ 1 ] [ <5T ] + 

ore*-.*- 3 -a^ - a..,) Pu lse[l] [T + 6T] + j«. M -4-« tt - 5 pulse [1] [2T + 6T1 + 

j«*-4. ««-,-«*.< > Pulse[l][3T + <5T] + J<-»-«-«»- 5 > Pulse[l][4T + <5T] + 
j<-n-.» Pulse[l][5T + <5T] + j(-«.-s * «»., > Pulse[l][6T + <5T] ) 

We select the imaginary 

0*c.N-5 ( j(a«.4*a N ., + a H . 3 .a M . 1+ a B ) ^^joj { M] + j(a M .4 + a N . 3 ♦ a N _ a ) pulse[0] [2 T + <5T] + 

J 0 "-* Pulse[0][4T + (5T] + Pulse[0][6T + <5T] + 

jti^.4*m-, --«-.) Pulse[l][3T + ST] + J^-«) Pulse[l][5T + 6T] ) 

ModulationValue[Pulse_] [{x0_, xl_, x2„ r x3_, x4_, x5_, x6_, x7_) ] 1ST 1 : = 
<!• % PulsefO] ISTI + (1. * i)-«>«4 Pulse[0] [2 T + 6T] + " 

(l. 1 I) x4 Pulse[0] [4 T + <5T] + (1. * I) ~ xS Pulse[0] [6T + <5T] + 
(1.* i)-^-x6-x7 E^isero] [8T+6T] + (i. » i)^*k3 + x4 Pulse[1] [T + 6T] + 
(1-* I) x3 Pulsetl] [3T + *5T] + (l.* I) ~ x6 Pulse[l] [5 T + <5T] 
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Moaulat±onValue[FiltPulse[8] ] [{-1, -1, -1, -1, -1, -1/ -1/ -1}][0.5T] // N 

-0.7540911 



ConvertToBitSeq[Depth_] [NumJ : = 
Module [{bit, n} , 
n = Num; 

ConvertNext : = (bit = Mod[n, 2]; n = Floor[n/2]; bit); 
( Table [ConvertNext, {i, 1, Depth}] //Reverse) /. {0 -1}]; 

ConvertToBitSeq[8] [2] 

{-1, -1, -1, -1, -1, -1, D -1} 

Table[i, {i, 1, 2 A 8>] // Map [Convert ToBitSeq [8] , #]&; 

Select [Table [i, {i, 1, 2 A 8)] // Map [Convert ToBitSeq [8] , #]&, (#[[5]] 



{ 



-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
1, 
1, 



-1, -1, -1, 1, -1, -1, -1}/ 

-1, -1, -1, 1, -1, -1, 1}, {-1, -1. -1, -1, 1, -1. 1, "D- 
-1, -1, -1, 1, -1, 1, 1}. {-1, -1, -1, -1. 1, 1* -1* -D' 



-1, 
-1, 
-1, 
-1, 
-1, 



-1, 

-1, 

-1, 

-1, 

-1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 



-l, -i; i, i,-i, i}, -i, -l, i, i, i, -i}, 

-i, -i, i, i, i, i), -i, -i, 1. 1, -i» -d -d. 
-i, i, i, -i, -i, i}, -i, i, -i/ 1. -i}- . , , , , 

-i, i, i, -i, i, i}, {-i, -i- -i, d -i* -d. - 1 ' x ' d x ' - 1 ' d, 

-l] 1, 1, 1, 1; -1), {-1, -1, -1, 1, 1, 1, 1, 1), {-1. -1* 1. -1. L - 1 ' - 1 ' 

-i; i, -i, i, -i, -i, i}, {-i, -i, i, -i, i, -i, i* -I), 

i -i i -i, i, i}' {-i, -i f i, -i, 1. 1, -i}, {-i. -i. -i. -i' d< 

1, -1, 1, 1, 1, -1}, {-1, -1, D -1. l.-l, 1* 1). (-1. -1- 1; 1; X ; "1; "D "D* 
1 1, 1, -1, -1, 1}, {-1, -1. 1. 1/1, -1. 1. "I)' ("I' X ' x ' X < 1 >' 

1, 1. 1, 1, -1, -1}, (-1. -1, 1, 1, 1. 1, -1, 1}, {-1, -1, 1, 1, I- 1' 1< , 
1, 1, 1, 1, 1, 1}, {-1, 1, -1, -1, 1, -1, "I, "Df (-1* 1* "I' -1' 1' - 1 ' "I' Df 
-1, -1, 1, -1, 1, -1}, {-1. 1, -1. -1- 1, -1, 1. 1>, nx 
-1, -1, 1, 1, -1, -1}, {-1, 1, -1, "I, 1/ 1, -1* 1)* {"!' 1' - 1 ' - 1 ' L 1 ' X ' " 1 >' 
-1 -l! 1, 1, 1. 1}, {-1. 1. -1- 1- 1. -1. "I- -D< ("I- !• 7 1 ' 1 ' lf ~ 1 ' V ' 

-1, 1, 1, -1, 1, -1}, {-1, 1. -1, 1, 1- -1, 1, 1}, {-!' If -1' 1' !' X ' - 1 ' 
-1, 1, 1, 1, -1, 1}, {-1, 1, -1, 1. 1, 1, 1, -1). {-1- 1' \* \' 1 \ 1 >- „ 

1,-1, 1, -1. -1, -1), {-1, 1. 1, -1, 1. "If "If !>' <-!• L X > D - 1 ' X ' ~ 1 >' 

1, -1, 1, -1, 1, 1}, {-1/1, 1, -1, 1, 1, -If -IK {"I- 1' If - 1 ' L 1* 1>* 
1, -1, 1, 1, 1, -1}, {-1, 1, 1, -1, 1, 1. 1, Df {-If If If 1; 1; -1; "I* -1>* 
1 1, 1, -1, -1, I-}, {-1. 1, 1. If If -If ^ "Df (-1. I', 1 ', 1 ' "V 1 ',J>' 
1, 1, 1, 1, 1, -1, -1), {-1, 1, 1. If If If -If Df {"I; LI' 1' 

1, 1, 1, 1, 1, 1, 1}, {If -If -If "If If -If -If -!>' {!' -1- - 1 ' - 1 ' - 1 ' X >' 
-1, -1, -1, 1, -.1. 1, -1}, {1. -If "If "If If "If 1' 1} ' ^ n _ - - - 

-1, -1, -1, 1, 1, -1, -1}, {1. -If -If "If If If "I- Df {!' - 1 ' 1' 1' 1.' "D' 

i, -i -i -i i, i, l. i}, {l. -if -if if if -if -if -Dftif r 1 ' TV'iVV 1 '^' 
i, -l, -l, i # i, -l, i, -l}, {i, -l, -if if if -if if I)- {1- -i. r 1 -, 1 ', 1 ', 1 ' *V," } ' 

i, -i, -i, 1, l, l, -l, l}, {i; -l. -1, if if if l'. -Df (if -l< -l-i' 1; i; D Df 

i, -i, i, -i, i, -i. -l, -d, {i. -i, i. -i, i, -1.1}' a* -v x - -V V "V i\ } ' 

1, -l. 1, -1, 1, -l, 1, Df {If -if If -if I- if -I- "I)' {D "V V V 

1, -1, 1, -1, 1, 1, 1, -1}. {If -If If -If If If If Df {!' -1; 1; D D "D T^' 

l, -l, l, l, l, -l, -1, D, {1/ -D if if if -if if -Df {!' -I* I* i' ' ~ V H^' 

X, -i, i, 1, 1, l, -l, -D, {l. -1, D D if if -D Df {D -D D l' D l, 1, -D, 

1, 

l, 
l, 
l ( 
l, 
l, 
1, 
l, 
1, 
l 
l 



-i; 1, 1, 1, 1, 1, D, {1/ If -If -If If "If -If "Df {If If "If "I; 1; "D -If Df 
i/-i -i, l, -l, l, -D, {l, l. -l, -I* l, -1, if Df {if I- -1' -I- I- D »V -i>' 
1, -i, -1, 1, 1, -1, 1), {1, 1. -1. -1, 1, 1/ D-D; {1; 1; ~ 1 ' ~D D D D D, 
1, -l, 1, 1, -1, -1, -D, {If If -If If If -D "If D- {If If -I- I- I- "D If -D' 
1, -l, 1, 1, -1, 1. l}, {l, 1, -l, 1, 1/ D -l, -D, {1, If -If If If I- ~ 1 ' Df 
1, -l, 1, 1, l, 1, -1}, {1. 1, - 1 » l f D if if Df {If I- if -I- 1 ' -l. "D^ -1} ; 
1, i, -i, 1, -1, -1, l}, {1/ 1/ l, -D 1, -1. 1, -D, {l.i.l' -if 1' -I. L I}, 
l, 1, -1, 1, 1, -1, -Df {if if if -1. D if -I- Df {D if if -l. Li. J * -D' 
1, l, -1, 1, 1, 1, D, {1, 1/ l, D 1/ -l, -l, -1}, {if If D if l. "D - 1 ' 1}' 
l, i, i, l, -l, 1, -l), {l, l, l, l, l. -l, l, D, {if l, l, 1, D D -I. - Df 
1, l, l, 1, l, -1, D, {1, D D 1. 1. if if -1}' {!' 1' 1' 1' 1. 1' D D) 
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xl ss Map[ModulationValue[FiltPulse[8] ] [#] [0.5 T]&, 

Select [Table [i, {i, 1, 2*8}] // Map [ConvertToBitSeq [B] , (#[[5]] — 1)&] ] 



{0.749266, 0.74922, 0.749173, 0.749219, 0.711482, 0.711529, 0.711482, 0 711435 
0.829796, 0.82975, 0.829703, 0.829749, 0.792012, 0.792059, 0.792012, 0 79196s' 
0.776885, 0.776839, 0.776791, 0.776838, 0.739101, 0.739147, 0.7391, 0.739054 ' 
0.802178, 0.802131, 0.802084, 0.80213, 0.764394, 0.76444, 0.764393, 0 764347 
0.759521, 0.759475, 0.759428, 0.759474, 0.721738, 0.721784, 0.721737, 0.72169 
0.819541, 0.819495, 0.819447, 0.819494, 0.781757, 0.781804, 0.781756, 0.7817l' 
0.78714, 0.787094, 0.787047, 0.787093, 0.749356, 0.749403, 0.749355 0 749309* 
0.791922, 0.791876, 0.791829, 0.791875, 0.754138, 0.754185, 0.754138, 6.754091 
0.749266, 0.74922, 0.749173, 0.749219, 0.711482, 0.711529, 0.711482, 0.711435 ' 
0.829796, 0.82975, 0.829703, 0.829749, 0.792012, 0.792059, 0.792012, 0.791965 
0.776885, 0.776839, 0.776791, 0.776838, 0.739101, 0.739147, 0.7391, 0 739054 ' 
0.802178, 0.802131, 0.802084, 0.80213, 0.764394, 0.76444, 0.764393, 0 764347' 
0.759521, 0.759475, 0.759428, 0.759474, 0.721738, 0.721784, 0.721737, 0.72169 
0.819541, 0.819495, 0.819447, 0.819494, 0.781757, 0.781804, 0.781756, 0 7817l' 
0.78714, 0.787094, 0.787047, 0.787093, 0.749356, 0.749403, 0.749355, 0.749309 # 
0.791922, 0.791876, 0.791829, 0.791875, 0.754138, 0.754185, 0.754138, 0.754091} 



Map[ll/2 - a/ 2 Erf [— ] U, xl] // Apply [Plus, #]& // #/128&; 



Ber[a_] : = Evaluate [*&4 8] 



Ber[0.01] 

0.495 



Erf [Infinity] 



1 



D[Erf [x] , x] 



2 E'* 2 
?Erf 



Erf[z] gives the error function erf(z). Erf[z0, zl] gives the generalized error 
function erf(zl) - erf(zO). 

Plot [Erf [x], {x, 0, 10}] 



10 



0.9999 
0.9998 
0.9997 
0.9996 



- Graphics - 



? FiltFulse 



Global 'FiltPulse 
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T := 3 / 812500 

Null 

?? T 

This is the symbol period 
?? T 

This is the symbol period 
T : = 3 / 812500 

Clear[T, xO, xl, x2, x3, x4 , x5, x6, x7] 

j(os-**aM^* «m-2 ♦ «H-1 ♦ «» > Pul se [0] [ 6T] + j<«« -*♦««-»♦*»- 3 > pul SO [ 0 ][ 2 T + 6T] + 

j**-* puisefO] [4 T + «5T] + j<- a »-*> Pulse[0] [6 T + <5T] . + 

j(-a H _ 5 -a M . c -a B _ 7 ) Pulse[0][8T + <5T] + j(a»-4 ♦ a M . s ♦ a« . i ) pulse [1] [T + <5T] + 
j(a..i) pulse [1] [3 T + <5T] + Pulse [1] [5 T + <5T] /. { a H -* xO, a M -i -» xl, 

a H . a x2, a N _ 3 x3, a N _ 4 -+ x4, a M _ 5 x5, a H „ 6 x6, a N - 7 x7) ^ 

(1. p u i S e [0] [<5T] + (1. i)«2+x3 + x4 p ulse [ 0 ] [2T + <5T] + 

(1. I) x4 Pulse [0] [4T + <5T] + (1. I)"* 5 Pulse [0] [6T + 5T] + (1. i)"* 5 '* 6 -* 7 Pulse [0] [8T + (5T] + 
{1. x) xl ***** A Pulse[l] [T + <5T] + (1. I) x3 Pulse[l] [3T + <5T] + (1. I)" x6 Pulsefl] [5 T + <5T] 

j{a K .j) pulse [1] f3 T + <5T] / . { a N -> xO, a s _i -+ xl, a N „ 2 -» x2, a N _ 3 x3, 
a N _ 4 x4, a N _ 5 *5, a„- 6 -+ x6, a N _ 7 -» x7 } 

(1. I) x3 Pulse [1] [3 T + <5T] 

(1. 11 I)* 3 Pulse[l] [3 T + 5T] 

(1. I) x3 Pulse[l] [3T + (5T] 
? J 

J = ° 



y this paper we study Laurents paper. 

™7new^d7a S that * " P ° SSible t0 C ° nStr " Ct PhaSC modu,ation *>V *e superposition of amplitude modulated pulses. 

We can separate the method of the superposition from the shape of the amplitude modulated pulses. We can then 
use different pulses of the same duration and /or fewer of them to construct modulation schemes that have desirable 
properties in the spectral and approximate to constant envelope schemes. 

In particular we can use pulses with narrower bandwidth (but of the same support) to generate a modulation scheme that 
tias a narrower bandwidth with a corresponding increase in the variation of the envelope of the signal . 

oz^JLJ^^ - JJL A /y^i^t Li 
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this paper we effectly study Laurents paper. 
We develop the modulating function that is used in GSM first to use as an example to demonstrate Laurent's idea. 
3 



812500 
BT : = 0.3 

V^og[2] 

o : = 



2 7TBT 

V 2 7T a T 



Ideally we would and did define the effect of the convolution of h[t] by the formula below. However, this version of 
Mathematica gives an error. / 



ct h[t - t] .. 
h f ii [t_] := Release [Module [{r}, J ^ - cStJJ 

We get round this by numerically fitting an interpolation function. 

h £11 [PulseWidth_] [t_] := Module [ {xl, x2, x3}, xl = 

rfh[tl-r] PulseWidtli P ulseWidth T . 
<N[#1, 40] t) [Table [{tl, t dr} f {fcl, , . ^)JJ 

T 

x2 = Xnterpolation[xl] ; x2 [t] ] 
Plot[h[tT], {t, -5, 5}, PlotRange -> All] 

- Graphics - 

In this section we define 0 L , t . the phase modulation function with a pulse width of length L at time t 

1 

Modulationlndex : = — 
2 

We use the definitions in the paper 

ff : = N [Modulationlndex 7r] 

C : = Cos [*] ; 
S := Sin [ft] ; 
j := ( e* * / / Chop) ; 
M := 2 1 *" 1 '; 

Clear [PhaseAngle] 

PhaseAngle [t_] /; t i 0 : = 0 

PhaseAngle [L_] [t_] /; tiLT:= » 
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") PhaseAngle [t_] : = 

Pha S8ta gle[L][tJ = Module[{xl, x2, x3, x4, x5, x6}, xl = h m [3 LT] [ tl - ±1], 

X2=Table » t2 ' *£ 3 T =valuate[x a] -tl}. {«, 0, * », _!_}], interpolation [x 2] [ 

RuleDelayed: :rhs : Pattern t_ appears on the right-hand side of 

rule 0 L _. t _ (^ it _ = Module [{xl, x2, x3, x4, x5 ( x6), xl = h m [3 LT] [tl - iilj, 
x2 = Table[ «1» ] ; lnterpolation[x2 ] [t] ] J . 

When running the notebook we fix the value of L here and the rest of the graphs and functions use this value . 
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L : = 8 ; 

phasepoints = Table [{t, PliasoFunctioix[t T] ) // N, {t, 0, I, , 1 / 40} ] , 
1,1 stPlot [phasepoints, Plot Joined -> True] 

NIntegrate: mcvb : NIntegrate failed to converge to prescribed accuracy after 7 

recursive bisections in tl near tl = -3.68196x10" . 
NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 

recursive bisections in tl near tl = -7.4914x10" . 
NIntegrate: mcvb : NIntegrate failed to converge to prescribed accuracy after 7 

recursive bisections in tl near tl = -5.43232x10" . 

G FCrther : output of NIntegrate: :ncvb will be suppressed during this calculation. 

NTnt e arate--slwcon • Numerical integration converging too slowly; suspect singularity, value 
NIntegrate .slwconwumer a integrand, or insufficient WorkingPrecision. 

°A your integrand is oscillatory try using the option Method->Oscillatory in NIntegrate. 

Nlntearate-'slwcon - Numerical integration converging too slowly; suspect singularity, value 
«f iAteoration i? 0, oscillatory integrand, or insufficient WorkingPrecision. 

• ?S your integrand is "oscillatory t£ using the option Method->Oscillatory in NIntegrate. 

WTn^mrate.-slwcon - Numerical integration converging too slowly; suspect singularity, value 
NIntegrate. .siwcon.wumer atory integrand, or insufficient WorkingPrecision. 

?£ you/inte^ -in^ the' option Method->Oscillatory in NIntegrate. 

G F^her : output of NIntegrate: = slwcon will be suppressed during this calculation. 
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1.25 
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0.75 
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2 4 
- Graphics - 

We can define an phase modulated signal to be 
S t = E 1 Z~-— a ° ^-° T 

As can be seen from the graph 0 t is 0 for most of the time and is nearly O after LT . We denote by * L , t the pulse that is 
set to these values 
i.e <t> hlt = 0 t<0; 

0l, t = 0t Otherwise 

S t = E I2 ~=— ffB * u, - nT 
We now consider t = N T + t where 0 < r < T 

S t = E IS "=— a * * L - T+(N_o)T 

Making use of the fact that 0 U , t = 0 t < 0 ; 
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) S t = E 1 ffn *L,r + <N-n)T 

S t = E 1 ^n=-« ff « *L.T + ( N-n)T £ ! N -L ♦ 1 °« *Ut ♦ (N-n)T 

Making use of <£>l , t = ^t^LT we get 
And using the notation for J : = 0 
Making the substitution n = N - i 

S t == J^-i^ E IX ^0 ,tt N-.^L,r + IT 

Express this as 



1 = 0 ' 



Since we have 



(Sin[A-B] Sin [B] ( e* x - ) \ 

I sl»[A] + e il^If Sin ^ - ^ L J / / Simplify 



E -IB 

And 
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) fin[A-E] + ^ Sin[B]_ ^ ^ // simplify 

[ Sin[A] Sin[A] 2j ) 

E IB 

We can make the following substitutions S := Sin[<J>] 

s = n 1 [ Sm[<P ~ 0l,t * it1 + r»- Sin[fc>T+iTl } 

1 1 = 0 Is S / 

We can make the substitution 

t_] /; 0 < t < LT := tf> L ,t 

t„] /;2LT > t i LT := 8 -^t-n 
We need to put this to avoid the function being extrapolated 

, t_] /; ! (0<t<LT) && 1 (2LT>t iLT) s = 0 
Clear [torn, xl] 

torn = <xl[t_] = (A [L # tT]; Table [{t, xl[t ]} // N, {t, 0, 2 1* , 1/40}]) 
ListPlot [torn, Plot Joined -> True] 

Now define 

Clear [ Laurent S] 

LaurentS[L_] [n_] [t_] := Sin[ tfr[L, t + n T] ] / S 
Plot [Laurent S [L] [2] [tT] , {t, -2, 10)] 

- Graphics - 

s, = J^«- n 1 f aB[ ^ lt * LTl + r-. ^Ilili ] 

i = o ^ S S / 

T N-L „ L-l f Sin[^t-NT + IT + LT] _„ N Sin[^ t - N T4iT] ^ 

S| = jZ2i.«. L ri\ LaurentS[i + L - NJW + J"""' LaurentS[i - N][t] ) 



i = 0 



We try to get an insight into V ( LaurentS[i + L - N][t] + J 0 *-* LaurentS[i - N][t] ) 
To this end we try defining 

Temp[L_, N_] [t_] : = fj (^[i + L-N] [t] + S[i-N] [t] ) 
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) (Temp[4, 0] [t] // Expand // Chop) /. Complex[0, 1.] -> Complex[0, 1] 

ia ~e?:r a :y a ° S[ ° ] [t] S[1] [t] S[2] [t J S W ft] + I 0 j -°- J * a -> 5(1] [t] S[2) [t] S[3) [t] ftl 
" " S 1 ^[°] ttj S[2) ft) 5(3] [t] SIS] It) *I--»*a. , 5[2] [t l j j 3 ^ i 4 [tl i 5 [tl i 

I- 5[0] [t] 5[5] [t] S[6] ft] 317] [t] + 5[4] [t] S[5) [t] i[6] [t] i[7] [i] ' ' 
In the next section we define the 

Notation[a Xi L_, K ^ii_ » AlphaKI[LL_] [K_, ii_]] 
Clear [AlphaKI] 

AlphaKI [LL__] [K_, i_] /; 0 < i < LL&& 0 ^ K < 2"- 1 : = 
Module [{xl, x2, x3, KNum), xl t= (x2 = Mod[KNuni, 2]; KNum = KWU1 ^" x2 ; x2 j . 
KNum = K; x3 = Table [xl, {i±, 0, LL - 1} ] ; x3 M ] 

Check the function. 

Sum[2*(ii - 1) AlphaKI [6] [13, ii] , {ii, 1, 6 -1}] 
13 

We can now define the LaurentC functions. 
Clear [LaurentC] 

» 

General: rspelll : Possible spelling error: new symbol name -LaurentC" is similar 
to existing symbol' "LaurentSV oa«u.j.« 

LaurentC [L_] [KJ [t_] /; 0 <: K < 2 L : = 
LaurentS[L] [0] [t] J"J LaurentS[L] [ii + L AlphaKI [L] [K, ii] ] [t] 

iiol 

LaurentLK gives the support of the C functions 

LaurentLK [L_] [K_] : = 

Module[{xl), xl = Table[L.(2 - AlphaKI [L] [K, ii] ) - ii, {ii, i, L - 1) ] ; Min[xl] ] , , 
Table [LaurentLK [L] [ii] , {ii, 0, 2 t - 1 -l}] 

{ !' I' ?' §' I' f' f' f' 4 ' 4 < 4 ' 4 < 4 - 4 ' 4 - 4 < 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 

2 2 2 2 2 2 I' V V V V V ?' ?' ?' ?' ?' ?' 2 ' 2 ' 2 ' 2 ' 2 ' 2 < 2 < 2 < 2 

1 1 ?' V ?' ?' ?' ?' ?' ?' J' }' lf lf X ' X ' X ' X ' 1 ' X ' X ' X ' 1. L 

i ?" i' i' n' X ' 1 ' X ' lf X ' X ' X ' X ' X ' 1/1* 1- 1. 1. 1, 1, 1, 

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, i, 1, l f l, l, i) 

CFunctionPlot[K_] := Module [ {xl, x2> , xl [t_] = LaurentC [L] [K] [t T] ; 

X L = ^ h , l6l } t: ' ^^H' t* 1 * °' LaurentLK [L] [K] , 1/40}]; ListPlot[x2, 
PlotJoined -> True, PlotLabel -> String Join [- \n \n C = », ToStringriill 
FormatType -> StandardForm, PlotRange -> All] ] 

Plot at most the first 7 C functions 



Table [CFunctionPlot[ii], {ii, 0, Min[6, 2 L " 1 -1]}] 
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0.0014 
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0.001 
0.0008 
0.0006 
0.0004 
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\ { - Graphics - , - Graphics - , - Graphics - , - Graphics - , - Graphics - , - Graphics 

- Graphics - } 

In this section we define the modulating function. We intend to build the modulator from first principles. 

V^-J H-l 
Z *K*0 



~~ 1 J* K * n C K [t -nT] 

J K*0 



Change the order of summation 

H-l ^ 

V J] O*^ C K [t-nT] 

Given that we are interested at t = N T + r 

M -1 ^ 

Change variable n -> n' + N 

M -1 m 

} ^ J***'*" C K [r -n' T] 
Notice that the support of C K is positive only. 

M -1 0 

y 1 2 j*K.n* ♦ n Cjc [r - n' T] 
Using LaurentLK as 



0 

^^n' o-LnurentLKfL] [Kj ♦ 1 



Changing the sign of n 

Laurent UC fL] [K] - 1 



2 



jA*, ».« c K [r + n T] 



Let the value of A 0 , o be specified. 

Then we can calculate the value of A 0 , -1 ,A 0 , - 2 > A o, -3 using the following algorithm 
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) BitSeq = Table [1, {i, l, 20}] 

{i, i, l. i, i, i, 1, i, 1, i, i, i, i, a , 1# 1# lf lf lt 1} 

ANKXnitialStateSetUp[I._] [K_] [XnitBitSeq_, AccumulatedPhase ] := 
InltbltS?: ^nSiSeq?' **' aCUPha8e ' ^^bitse,), " 
acuphase = AccumulatedPhase, 
UpdateSeq : = 

Module[{}, xl = acuphase - Sum[initbitseq[ [i ]] AlphaKiri.] [K, il ti ± x. ni 

TablefunL; Q aCUP ?. aS « " ^rsttinitbitsa,]; initbits^q = Rest [ inltLi'tseq? A\ 

Table[ UpdateSeq, { i# 1, LaurentLK[L] [K] } ] ] "oitseqj, xxj 

With the BitSeq we have we expect a decreasing sequence starting at 0. 
ANKInitialStateSetUp[L] [0] [BitSeq, 0] 
{0, -1, -2, -3, -4, -5, -6, -7, -8} 

Testlnit = Table [1, {i, 1, L}] 
{1, 1, 1, 1, 1, 1, l, 1} 

AKN[L_] [K_] [{State., AccumulatedPhase J ] : = 
AccumulatedPhase - Sum [State [ [i + 1] ] AlphaKI [L] [K, i] , { i, 1, L _ 1} j 

Options [Modulator] : = {StartingQuadrant -> 0, InitialState - Table U fi 1 20n 
Sainplinglnterval - T/32, Nu^nberOf Curves - 4, ModulatingPuIse^ LaurektC } ^ } ] ' 

BitSeq = Table[l, {i, 1, 20}] 

{1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, !, lf lf lt lt lf 1} 

Modulator [L__] [Bit Seq_, Opts ] : = 

Module [ 

state = InitialState /. {Opts} /. Options [Modulator] , 
x3 = StartingQuadrant / . {Opts} /. Options [Modulator] , 
X4 = SamplingXnterval /. {Opts}/. Options [Modulator] ; 
Pulse = ModulatingPulse/. {Opts} /. Options [Modulator] , 
Curves = (Numberof Curves /. {Opts} /. Options [Modulator] ) - 1, 
AccumulatedPhase = x3 - BitSeq[[l]], J ; ' 

seq = BitSeq; 
Table[AKNState[K] = 

x5 .^odul«?/ 1 f ta, ; e t et0P[ T I ' ] . [K r ] [ ! tate ' *«umulatedPhase], {K, 0, Curves}], 
x5 Module[{}, state = Join [ {First [seq] } , Drop [state, -111, 
AccumulatedPhase = AccumulatedPhase - First rseal • 
saq = Rest [seq], 1 ,J ' 

Table [AKNState[K] = Join[ {AKN [L] [K] [ {state, AccumulatedPhase} ]} , 
Drop[AKNState[K], -1]], {K, 0, Curves}],- 1U ' 
*6[r_] = sum[sum[ (J) tli.ijj pul-e + ± 

n-M r JV, °' 1 . LaU « entI * K[r ' ] CK] _1 >1' < K ' °' Curves}], 
Table[x6[r], {r, 0, T - x4, x4}]], 1 

Table[x5, {He, 1, Length [BitSeq] } ] //Flatten] 

General: : spel 11 : 

Possible spelling error: new symbol name "state- is similar to existing symbol -State". 

^"simiU?^ 111 • f? ssible selling error: new symbol name •'AccumalatedPhase- is 
similar to existing symbol "AccumulatedPhase- . »^«araase is . 



La u ren tPaperPaten ts.nb 



1 1 



RandomBitSeq = Table [Random [Integer, {0, 1}], {i, 1, 40}] //Map[# (-2) + 1ft, #]& 
{1, 1, l f 1, 1, 1, 1, -1, 1, -1. 1, -1; 1- 1- "I' ^ ~\ - 1 ' - 1 ' - 1 ' 1 ' 

i ( i, l, a, -l, -i, -i* - 1 * !< - 1 * x > " 1} 

Modulator [L] [BitSeq, Samplinglnterval -» T/10]; 

Modulator [I*] [BitSeq, Samplinglnterval -+ T/10, KumberOf Curves -> 8] ; 
Modulator [L] [BitSeq, Samplinglnterval T/10, NumberOf Curves -> 1] s 
torn = Modulator [L] [RandomBitSeq, Samplinglnterval T/10, NumberOf Curves -> 2] 
ListPlot [Im[tom] , Plot Joined -> True] 




400 



- Graphics - 
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) ListPlot[Abs[tom] , PlotJoined -> True] 

1.002 



1.001 



0.999 
0.998 
- Graphics - 



M 



W 



M 



M 



3 DC 



400 



{Re [torn], Im[tom]> // Transpose // Lis tPlot [#, PlotJoined -> True, AspectRatio 
-4- 




Graphics 



torn = Modulator [L] [RandomBitSeq, Samplinglnterval - T/10, NumberOf Curves -> 1 
ListPlot [Im[tom] , PlotJoined -> True] 




400 



Graphics 
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ListPlot [Re [torn] , Plot Joined -> True] 




400 



- Graphics - 

ListPlot[Abs[tom] , Plot Joined -> True] 



1.06 
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0.98 
0.96 
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4 00 
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- Graphics - 
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{Re [tea], Im[tom]> // Transpose // Li.tPlot [#, PlotJoined -> True, AapectRatlo -> i ]& 
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- Graphics - 





Calculating the bandwidth of the signal 

xll = Table [Laurent C [L] [0] [t T] , {t, 0, L + 1 - 1/200, 1/200}]; 
((Drop[xll, -1] -Drop[xll, 1]) 200 ) *2 / 200 // Apply [Plus, *]* 

1.29684 

Here the pulse 6T is just stretched to 8T to see the performance 
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PhaseFunction[t__ ] = 4>s,*t 

phasepoints = Table[{t, PhaseFunctionft T] } // N, {t, 0, 6, 1/40}]; 
ListPlot [phasepoints. Plot Joined -> True] 

NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 1 . 91796 x 10 -7 . 

NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = -5.52869x10" . 

NIntegrate- :slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision 
If your integrand is oscillatory try using the option Method- >Oscillatory in NIntegrate. 

NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = -1.96696x10" . 

General :: stop : , . . , . . 

Further output of NIntegrate: :ncvb will be suppressed during this calculation. 

NIntegrate- :slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision. 
If your integrand is oscillatory try using the option Method- >Oscillatory in NIntegrate. 




atory in NIntegrate. 

^Srther" output of NIntegrate: : slwcon will be suppressed during this calculatic 



ion. 



- Graphics - 

TestPulse[8] [0] [t_] = LaurentC [6] [0] [t 7 / 9 ] 

l.S±n[*[6. ^ikoO^^ Sin[*[6,- — *— + 2*.)] Sin[*[6. ^00^)) 

^l*l*> ToAlT^ll 31 ^^ 6 ' T6B00 + ^]l sin W 6 ' ^ 

Plot [TestPulse[8] [0] [t] , {t, 0, 9T} # PlotRange -> All] 



- Graphics - 

RandoxoBitSeq = Table [Random [Integer, {0, 1)], {i, 1, 100}] //Map[# (-2) + 1&, #]& 

{i, i, -l, -l, -l, 1. 1. -l. i. -i, i- -i. i. - 1 ' - 1 ' 1- ; lt r\ t "V lf 1# lf 
i. -i, i, -i, i, i, i, -i. -i. -i. i, -i. v V V V ~n V V 'i i " i i ii 

1, 1, 1, -1, 1, -1, -1, 1,. 1, 1, -1- -1* 1* 1. !' !< X ' X < 1 ' X ' _1 ' X ' * 1 ' X ' " 1} 

torn = Modulator [L] [RandomBitSeq, Samplinglnterval -> T/10, NumberOf Curves -> 1, 
ModulatingPulse TestPulse] ; 

ListPlot [Im[tom] , PlotJoined -> True] 



* - Graphics - 

ListPlot [Re [torn] , PlotJoined -> True] 



- Graphics - 
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) ListPlot [Aba [torn], PlotJoined -> True] 

- Graphics - 

{Re[tom], ln ,[tom]} // Transpose // ListPlot [#, PlotJoined -> True, AspectRatio -> 1], 

- Graphics - 

In this section we consider filtering the pulse 
Needs [ "bessel * n ] 

fil = BesselFilter[8][LowPass3dB[l/T 2 Pi 0 . 3] ] [Fre«juencyResponse] [s] ; 

Plot [20 Log [10, Evaluate [ (fil /. s -+ 2 Pi I f / t) // shall re n r>\ 
GridLines -> Automatic, AxesLabel -> { -1/T-, -dB«T] 



- Graphics 



Res[t_] = InverseLaplaceTrans form [fil, s, t] 

2027025 ((9.36812 + 0. I) E" 897172 * fc Cos[139301 t] - 

(13.7944 + 0. I) E -83567 2 .t cos[420044. t] + (4.73111 + 0. I) e -701358, fc Cos f 708768 M 
(0.304883 + 0. I) E-«»»-t cos [1. 02016x10' t] + Cos[708768. t] - 

(29.1264 + 0. I) e- 897172 -' Sin[139301. t] - 

(9.89398 + 0. I) e~ 835672 * fc Sin[420044. t] - (0.321205 + 0. I) E^ss. t sin , 708?68 
(0.375154 + 0. I) E-^ssae.t sin [1 . 02016 x 10 fi t] ) S«i[708768. t] 

Plot [Res [t T] , {t, 0, 6}, PlotRange -> All] 
- Graphics - ■ 

filTaps = Table [Res [t T] , {t, 0, 6, 1/20}] // Chop; 

We estimate the group delay to be 1.8T 

Length[filTaps] 

121 

pulse =Table[LaurentC[8][0)[tT], {t, 0, 9, 1/20}]; 

General: ispelll : 

Possible spelling error: new symbol „a,e "pulse" is similar to existing symbol "Pulse". 
pulse2 = Table [LaurentC[8] [1] [tT], {t, 0, 7, 1/20}]; 
Needs [ " Convolve » " ] 



LaurentPaperPatents.nb 



FiltoredPulee = T/20 MapConvolve [pulse, filTaps] //.Take[#, {34, Length[#] )]& 

{1 73731x10"" 5 35887x10"", 1 . 60413 x 10" 14 . 4 . 66278 x 10' 14 , 1 . 31688 x 10" 13 , 
3*61578xl0- 13 ! 9. 65729X10" 13 , 2 . 51044 x 10" 12 . 6 . 35507 x 10" 12 , 1.56747x10"", 
S^eeegxlO" 11 , 8.83867 x10-", 2.02275xl0- 10 , 4.51958xl0- 10 , 9.86455 xl 0 - , 
2 'l0426xl0- 9 , 4.38922X 10" 9 , 8 . 95697 x-10" 9 , 1 . 78914 x 10" 8 . 3.4999xl 0 - 8 6 . 70842 x 10" 8 . 
1 26054xl0" 7 , 2.32319 xl0- 7 , 4 . 20166 x lO' 7 , 7.46062x10"'', 1.30124x10" , 
2'23O35xl0" 6 3 75853 xl0" 6 , 6.23xl0- 6 , 0.0000101618, 0.0000163174, 0.0000258049, 
n'nnon!o2068' 0 000061746, 0.000093496, 0.00013964, 0.000205786, 0.000299337, 
S'000429924 0 6o0609895, 0 000854853, 0.00118424, 0.00162196, 0.00219697. 

0 0029439 0 00390362, 0.00512368. 0.00665873, 0.00857068, 0.0109288 0 0138093, 
n'n?75^2' o 021476 0 0264451, 0.0323001, 0.0391407, 0.0470666, 0.0561754, 
S"o665603' 0'0783079 6 0914948. 0.106186 0.12243, 0.140261. 0.159691 0 180711, 
n' 20329 0 22737 0 25287. 0.279682, 0.307673, 0.336686, 0.366539, 0.397029, 

0 459013 0 490013, 0.520668. 0.550708, 0.579856, 0.607837, 0.634379, 
S*M922 0 682107 0 702805, 0 721098. 0.73679, 0.749713, 0.759727, 0.76672, 
2"t70«16 6 771369 6 768969, 6.763442, 0.754846, 0.743274, 0.728852, 0.711736, 

1 sllllV C \ eillll: 0 elllll: 0.62035. 0.592964, 0.56429, 0 . 534607 0 -504193, 

0 442277 0 411305, 0.380656, 0.350558, 0.321221, 0.29283, 0.26555, 
S'SJII?!' O'itllll: Si 191636 O. 169 9 36 0.149791, 0.13122, O.II4217 00987613. 
00848118 6 0723143, 0.0612016, 0.0513967. 0.0428146. 0.035365, 0.0289542, 
n n5?ABT 0 0188687 0 0150065 0.0118109, 0.00919671, 0.00708412, 0.0053992, 
S' 00407439 °'0 487?', 0 TooIleSll. 0.0016848. 0.00125761, 0.000951417 0.000736739. 

nnno^R9309 6 000489543 0.000422011, 0.000374899. 0.000339478, 0.0003096, 
0"000281223' o'o00251974 0 000220764 0.000187441, 0.000152497, , 0.000116824. 
oononRl^lis 6 0000477044 0.0000164566, -0.0000113198, -0.0000349206, 

0 0000538729 lO 0000679391. -0 0000771062, -0.0000815622, -0.0000816651, 
:2"oSSo77?073; -0".oS007oi77, -0 000061222, -0.0000496142, -0.0000367186, 

n'nnnn27ififiR 9 53584xl0" 6 , 3 . 66817 x 10" 6 , 0.0000160205, 0.0000271838, 

0 oSSSailol! 0 0600450345. 0 0000514769, 0 0000562268. 0.0000593363, 0.000060909 

S'SSSSI"?!!: o^oweoolsl 0:0000579702: o.oooossosag. o ooposisoTe. o. 0000475036, 

00000432203 0 000038813, 0.0000344153, 0.0000301387, 0.0000260716, 
00000222801 0 ! 0000188096, 0.000015687, 0.0000129228, 0.000010514, 
8'4469xl0- 6 , 6.69967xl0" 6 , 5 . 24481 x 10" s , 4 . 05145 x 10" 6 , 3.08723x10" , 

2 3199xl0" 6 1 71859xl0" 6 . 1 . 25466 x 10" 6 , 9. 0236xl0" 7 . 6.39109x10 , 
4'45602xl0" 7 , 3.05722xl0" 7 , 2 . 06318 x 10" 7 , 1 . 36898 x 10" 7 , 8.92721x10" , 
5'71866xl0" e . 3.59686xl0" 8 , 2 . 22017 xlO' 8 , 1 . 34418 x 10" 8 , 7.97812 x10" 

4 . 63953 x 10" 9 , 2 . 64191 x 10" 9 . 1 . 4722 x 10" 9 8 . 02333 x 10"» 4 . 27381 » x 10" 

2.2234x10"", 1.12886x10"". 5.59276x10-". 2.70717 x10"", 1.28211x10 .. 

5.92683 x10"", 2 . 65259 x 10" 12 , 1.. 14156 x 10"» . 4.73562x10" , 1 . 9432.9 x 10" 

7 95416X10" 14 , 3.14517X10" 14 , 1 . 23089 x 10" 14 , 4.73302x 10"". 1.66595x10 

4 65283x10-". 1.93314x10"", 9 . 87738 x 10" 17 . 2.31472x10"". - 1 . 50823 x 10" 1 

-1 01228X10" 17 , -6.49877x10"", -2.08455x10"", 1.80177x10"". 1.70062x10" , 

1 9917x10"", -1.14008x10-". - 1 . 66045 x 10" 20 , 3 . 7877 x 10"" , 1.47752 x10" . 
l'94753xl0- 21 , 3.72298 x 10" 22 . 2 . 0 2 5 3 5 x 10" 22 , 1 . 36007 x 10" 22 , 6.72514x10" 

7 '. 07296 x 10"" , 6 . 84862 x 10" 23 . 3 . 6 8 9 1 6 x 10" 23 , 2 . 6013 x 10" 24 . - 1 . 51251 x 10" 

-5 51867x10"", 7.52915X10" 26 . 5 . 8 1 5 8 4 x 10" 25 , 1 . 75371 x 10" 25 , 1.22392 x 10 , 0} 



LauretuPaperPatents.nb 



18 



FilteredPulse. . t/ 20 HapConvolve [pulse*, filTaps] //Take [*, {50 , Leilgth[#] 

^.•sii":^ 7 ; f- o8o6xi °^ ™*io-. i.^xo-., 

0.0000313562 0;0000456 57 2/0'0^ 0.0000212804 
0.000251616, 0.000342391, 0 000460925 c T 'nnnc??2^ 48 n 12 ^° " 000131493 ' 0.000182902, 
0-00136282, 0.00174186 0 0022o2l« n 'nSoS2?fi 39 « 32 ' °-°00809185, 0.00105553/ 
0.00618453: 0.007391o!,' S.'SSIIm"' " 3 ' 0 °n ?^^° 5 o 0 00421186, 0.00512763. 

0.0178876, 0.0201213 0 022«ng n n^Sl? « ;,J; 0119544 ' 0.0137895, 0.0157714 
0.0345256 0.0367857,' S.'SllloBs' 0*04085^ °n°UllZh °: 02974 19. 0.0321652, ' 
0.0463073, 0.0469636 0 0473081 o'oJ?S!«' n'n^n?«^' S^ 107 ' °-°453496, 
0.0443919, 0.0429671 0.0413131 0 0394S9 'n n^I?^ 3 'n°A 046452 5 ' °-0455611. 
0.0306913, 0.0283283 0 0259611 o'o2? 45 «i A°n? 4 |f 4, °-°3S2785, 0.0330189, 
0.0149823, 0.0130987 oioilffoi' o'ooII^r' S " nnpo 2 ,^ Q ° ' 019112 ' 0-0169912, 
0.00476282, 0.00386365, 0.0030899 0 lltll^ ' Ton?*^ o° n " °°^ 6697 ' °-00579S18, 

0. 000765388, 0.000536896. 0 000363597 T o 23^ n^\^^°^ 142752 ' 0.00105885, 

: p!ii:J:,TOiJ!;S:j^ 

1. : SSJSI; V^F •• ! !: ™": :SS °"-1<: ::«388: 

-0.0000174907, -O.OOOO^Sl Jo^M^ -0.0000134543, 



3.57494xl0* 6 
7.24522xl0" 6 
5.41573 xl0" 6 
2.10471xl0" 6 , 
4.82884 xl0" 7 , 
9.65984 xl0" 8 , 
8.74155xl0 -9 , 
7.73129xl0~ 10 
4.40884 x 10" 11 
1.6239 xl0" 12 , 



4.95319xl0 -6 
7.16306x 10" 6 
4. 81423 X 1CT 6 
1.70583 x 10~ 6 , 
3.59854xl0" 7 i 
6.71324xl0" 8 , 
5.55593 xl0" 9 



5.98756xl0" 6 
6.89804xl0" 6 , 
4.20863 x10^ 
1.36297xl0~ 6 , 
2.64486xl0" 7 



4.54101X10' 10 , 2.62031x10 



4.59913 xl0- 8 , 3.10539xl0" 8 
3.47564xl0~ 9 , 2.13956 xl0" 9 



■10 



-2.1065xl0' 7 , 1.84943 xl0" 6 , 
6.69487xl0' 6 , 7 . 10243 x 10~ 6 , 
6.49223 xl0~ 6 # 5 . 98586 x 10" 6 , 
3.6208xl0" 6 , 3.06728xl0" 6 , 2 . 55958 x 10" 6 
1.07376X 10- 6 , 8.34162xl0- 7 , 6 . 3908 x 10' 7 ,' 
1.37044 x 10" 7 , 
2.06612xl0" 8 , 
1.296l3xl0" 9 



1.91715xl0- 7 ( 



1.3542 x 10* 8 , 



^ -~ , 1.47888X 10" 10 , 8.15113 xlO' 11 

2. 37341x10-", 1.26865x10-", 6.58958x10"" 3.34992x12-' 

. 7 -°_i 03 a 7 c, X , 1 1 ° ^- 33 °. 365<10 - 13 ' '-75051x10-, '-3. 22249x10-; 



-8.12875x10-", -7.89371X10" 14 



4.12557X10-, L9906-X1-0-; 1 . Wo68 x 10- 10 „ 4 ^Ul^^lU 

s "til: 0 -,i^ i;^ iii^iv;?- j- 92164 * 10 -;:- «.3i9«xi 0 -». 

3. 90929x10-, 1 S "Jo"-. 0} 2 X ' 7 " 33335 * 10 "" < 1-13717x10-, 

FiltPulse[8] [0] = { (Table [tT, {t, 0, 2 0, 1/20)1 // *n 1 1 

// Transpose // 

InterpolatingFunction[{{0, 0 . 00004929 2 3} } , <>] 

riltPulse[8][l] = {(Table [tT, (t, 0, 20, 1/20}] // N ) // 

e2> // Transpose // 

InterpolatingFunction[{{0, 0.000038953 8}}, <>] 



14 



-8.21475xl0" 15 
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Plot[FiltPulse[8] [0] [t T] , {t, 0, 9}, PlotRange -> All] 



0.6 



0.4 



0.2 




2 4 6 8 

- Graphics - 

Plot[LaurentC[8] [0] [t T] , {t, 0, 9), PlotRange -> All]- 




. - Graphics - 
Show[%, %%] 




Graphics - 
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Plot[FiltPul S e[8]tl]ttT], {t, 0, 7>..PlotRange -> All] 




1 2 3 4 5 6 7 
- Graphics - 

Plot[LaurentC[8][l][tT], {t, 0, 7), PlotRange -> All] 




- Graphics - 

Bandwidth [TestPulse_, Supp_] •= 

Bandwidth [FiltPulse [8] [0], 9] 

0.711659 

Bandwidth [LaurentC [8] [0] , 9] 

1.29684 

Bandwidth [FiltPulse [8] [1],*7] 

0.00338446 

torn = Modulator [L] [Table [1, {i, 1, 100}], Samplinglnterval T / 10 
NumberOf Curves -> 2, ModulatingPulse -> Fil^PuTs^ [? erVal T / 

^^^^U^S^S^ ^^nterval - T /10 , ^erOf Curves 



ListPlot[lm[tom], PlotJoined -> True] 



Graphics 
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ListPlot [Re [torn] , Plot Joined -> True) 

- Graphics - 

ListPlot [Abs [torn] , PlotJoined -> True] 

- Graphics - 

{Re[tom], Hn[tom]} // Transpose // ListPlot [#, PlotJoined -> True, AspectRatio -> 1]& 




- Graphics - 

torn = Modulatory) [RandoxnBitSeq, Samplinglnterval -> T/10, NumberOf Curves -> 1, 
ModulatingPulse -> FiltPulse] ; 

ListPlot [Im [torn] , PlotJoined -> True] 




Graphics 
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) ListPlot[Re[tom] , Plot Joined -> True] 




l/wv 



)0 



600 



l/WWWW 




000 



- Graphics - 

ListPlot[Abs[tom], Plot Joined -> True] 




- Graphics - 

{Re[tom], imitom]} // Transpose // ListPlot [#, PloMoined -> True. AspectRatio -> l ]& 




- Graphics - 

xll = TabletTestPulsetfi] [0] [tT], {t, 0, 7 - 1/200, 1/200}],- 
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(<Drop[xll, -1] -Drop[xll, X]) 200 ) *2 / 200 // Apply [Plus, #]& 

0.92632 

TestPulse[6] [0] [t_] = LaurentC[4] [0] [t 5/7 ] 

l.Sin[*[4, ^ w + -^]]sin[^[4. * **] ] Sin[*[«. gT^oO ^11 

Si»[*[4. if]] 



L := 5; 
PhaseFunction[t_] = fa.t? , 
phasepoints = Table [{t, PhaseFunction[t T] } // N, {t, 0, L , 1/40}]; 
ListPlot [phasepoints, PlotJoined -> True] 

NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 1 . 48053 x 10 -6 ♦ 

NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 1 . 12255 x 10" 6 . 

NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = -1 . 08808 x 10' 6 . 

further * output of NIntegrate: :ncvb will be suppressed during this calculation. 

NIntegrate: :slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand/ or insufficient WorkingPrecision 
If your integrand is oscillatory try using the option Method- Oscillatory in NIntegrate. 




1 2 3 4 5 

- Graphics - 
L : = 6 

TestPulse2[6] [0] [t_] = LaurentC[4] [0] [t 6/7 ] 

l.Sin[*[4. ^ I ^ F + ^]]sin[^[4, + ] Sin[*[4, ^00^)} 

Sin[#[4. " ]] 



torn = Modulator [L] [RandomBitSeq, Samplinglnterval - T/40, NumberOf Curves -> 1, 
ModulatingPulse FiltPulse] ; 

Save["ModulatorData.xn", {T, L, RandomBitSeq, FiltPulse, torn}] 
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ListPlot [Im[tom] , PlotJoined -> True] 




000 



- Graphics - 

ListPlot [Re [torn] , PlotJoined -> True] 

A/Wi 




L0 



lfWV 



00 



WWWW 




000 



- Graphics - 

ListPlot [Abs [torn] , PlotJoined -> True] 



4f 00 




Graphics - 



't 
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- Graphics - 

xll = Table[TestPulse2[6][0][tT], {t, 0, 7 - 1/200, 1/200}]; 



0 
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Needs ["Lauren tFunct ions x "] 

RuleDelayedrrrhs : Pattern t_ appears on the right-hand side of rule 

PnaseAnglelL_][t_] ( PhaseAng.le [L ) [t_] = Modulef {xl, x2 x3 , x4 x5 x6 } , «1» J ) . 

Keeds [ "LaurentNotationTest *" ] 

Needs: :nocont : Context LaurentNotationTest' was not created when Needs was evaauated. 

Information on the functions used can be obtained using help. 
Naities [ "Laurent Functions x *" ] 

{AXN, AlphaKI, ANKInit ialStateSe tUp, BT, FiltPulse h hFiltered miH.Kh^ t 



T := • 

812500 
BT : = 0.3 



Modulationlndex : = — 
2 

« ModulatorData.m; 
RandomBitSeq 

i: ^^^:,rA: i: ;S;V-;^- 1 ;V.V.V.V.V. fc'^c^. ,. .,.> 

M-l u _, 



tom2 = Modulator [L] [RandomBitSeq, NumberOf Curves -> 2, ModulatingPulse ^ FiltPulse] 
Li S tPlot[{Re[tom2], Xm[tom2]> //Transpose, PlotJoined -> 'True] 




- Graphics 



RandomBitSeq 
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Receiver [L] [tom2, StartingQuadrant -> 4] 
We write all the possible terms 

jAo.n-5 ( j(a»-4* <*n-s+ on -a ♦ on-i ♦ a « > Pulse [0] [ <5T] + 

J (o s -4*a». J .o„. a * «h-i) Pulse [0] [T + 6T] + ««-■> Pulse[0] [2 T + <5T] + 

Pu lse[0] [3T + <5T] + J*"'* Pulse[0] [4 T + <5T] + Pulse [0] [5 T + <5T] + 
j(-«»-s> Pulse[0] [6T + <5T] + j^-'- a -'> Pulse[0] [7 T + <5T] + 

j( -« H . s -a H . < -a.s-7) p u i SQ [0][8T + <5T] + jch. — N- a *--i ••h -h-i> Pulse [1] [ <5T] + 
J( ..-4.*.-.-«»-a*-»-i Pulse [1][T + 6T] + jO-4 > Pulse[l][2T + <5T] + 

JI...4...-.— -4i Pulse [1] [3 T + 6T] + Ji<-« — Pulse [1] [4 T + 6T] + 
j(-h-.> Pulse[l][5T + <5T] + j<-«-*«--'> Pulse[l][6T + <5T] ) 

/ j(a K -4 ♦ «h + °h-3 ♦ a«- l ♦ a w ) pulse [ 0 ] [ 6T] + 

j^.^ano-H.^a-.i) Pu lse[0][T + (5T) ■+ j(«,,.4-..o ♦ ««-») pulse[0][2T + 6T] + 
j(««-4-a*->) p u lse[0][3T + <5T] + J* M -* Pulse[0][4T + <5T] + Pulse[0][5T + <5T] + 
j(-««-s) PulseTO] [6T + <5T] + j(-«h- S -o*-«) P ulse[0] [7 T + <5T] + 

.,(-«...-«... —-i) P ulse[03[8T + ST] + Pulse[l][6T] + 

j(a M .4 ♦ on. j ♦ a N . 2 ♦ a N _ i - a N - a ) p U lse[l] [T + <5T] + 
j<a K _4 ♦ on -a ♦ + a« - a ) pulse [1] [2 T + 6T] + 

j(a N „ 4 *a N . s -cih-4) pulse [1] [3 T + <5T] + j(«»-4 - ««-s ) p u lse[l] [4 T + <5T] + 
jt">»-<) Pulse[l][5T + ST] 4- j(-OH.5 + a N _ 7 ) p u lse[l][6T + <5T] ) /. { a» - x0, 

j X o+xi+x2+x3*x4 Pu i se [0] [<5T] + j*i+ x2+x3+x4 Pulse [0] [T + 6T] + 
jx2 + x3 + x4 Pulse [0] [2T + c5T] + J x3+x4 Pulse [0] [3 T + <5T] + J x4 Pulse [0] [4 T + <5T] + 
PulserO] [5T + 6T] + J" x5 PulsefO] [6T + < 5T] + J-**-* 6 Pulse[0] [7 T + <5T] + 
J-X5-X6-X7 Pulset0 ] [8T + <5T] +J *o-x 2 .x3*x4 Puls e[l][6T] + Pulse [ 1 ] [T + <5TJ + 

jx2 + 2x3*x* p u i S e[lJ [2T + <5T] + J x3 Pulse[l) [3T + <5T] + J x4 " x5 Pulse[l] [4T + <5T] + 
J" x6 Pulse[l] [5T + 6T] + J" x5 * x7 Pulse[l) [6 T + <5T] 

Clear [XbsoluteValue] 

AbsoluteValue[Pulse_] [{x0_, xl_, x2_, x3_, x4„, x5_, x6__, x7_) ] [<5T_] : = 
Module [{al, a2, a3 , a4} # 
al - j* 0+ « 1+ ^* x3+x4 Pulse [0] [<5T] + 

j^x2^»4 Pu lse[0] [T + 6T] + a« a ** s «* Pulse[0] [2T+ <5T] + J x3+x4 Pulse[0] [3T + 6T] + 
j 314 Pulse [0] [4T + 6T] + Pulse [0] [5X + <5T] + J' x5 Pulse [0] [6T + M] + 

p ul 3 e[ o] [7 T + aT] + ^ Pulse[0] [8 T + ST] + j-*—**-** Pulse[l] [6T] + 

^i + x 3+ x4 puxsetl] [T + 5T] + jxa +2 x3 + x4 Pulse[ l] [2T + 5T] + J** Pulse [1] [3T+6T] + 
J**-* 5 Pulse Til r4T + i5T] + J"" 6 Pulse [1] [5 T + ST] + j'* 5 ** 7 Pulse [1] [6T + 5T]; 
a2 = ( Re[al]*2 + Im [al] - 2 / / ComplexExpand) /. {Im[x_] :> 0, Re [x_J => x> ; 
a2 /. x_ ? Number Q y_ -> Round [x] y] 

XbsoluteValuelnterf earingregions [Pulse.] [ 

{x0_, xl_, x2_, x3_, x4_, x5_, x6_, x7_} ] [<5T_] : = 
Module [{al, a2, a3, a4} r 
al _ J *o*xi*xa*x3*x4 Pulse [0] [fiT] + J x2 * x3+x4 Pulse [0] [2 T + ST] + J" x5 Pulse [0] [6T + 6T] + 
jlxs-xe-x? j^isero] [8 T + <5T] + J* 1 ^- 4 Pulse [1] [T + <5T] + J" x6 Pulse [1] [5 T + 6T] ; 
a2 = { Re[al]*2 + Im[ai]*2 //ComplexExpand) /. {Im[x_] :> 0, Re[xJ :> x) ; 
a2 /. x_? Number Q y„ -> Round [x] y] 
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AbsoluteValue[Pttls.[a] ][{-!, -1, -1, -1,-1, -i, _ 1# -l) H o.5T] 

x0[4] 2 + x0[12] 2 -2x0[4) x0[20] + xO [ 20] 3 - 2 xb [ 12 ] xO[28] + xO [ 28 ] * + 2 xO f 4 1 x0r 3K l 
X0J20] xO[36, + x0[36]> + 2x 0[ 12] xO [44 ] - 2 xO [ 28 J x0[44] + xO [ 44 ] » xO [i Txi f 5 
2 x0 20,xO[S2]-2x 0[ 36, xO [ 52 ] + xO [ 52 ] 3 - 2 xO f 12 xO 60 + 2 xO [ 2 8] xO Teo , 1 

xO 4, x0[60] + x0[60] + 2xO[4,xO[68J -2x0(20] xO [68] U xO [L Lo^T- 
2x0[52] X0[68] + xO [ 68 J * + 2 xO ( 12 ] xl[4] -2x0[28] xl[4] + 2x0[44] xl[4] 
2x0[60]xl [ 4 ]+ xl[4 ] 3 - 2 x0 [ 4]xl[12] + 2x0[2(i ) xl [ lV-2x0[i6]xlU2 + 
2 Cull X ][ll\- 2 *°W ^[12] + x 1( 12] 2 + 2x0[12] xl 20] - 2 X0 28JX 20] + 

xl 28]^2x0[12] xl[36] -2xO[28] xl [ 36] ♦ 2 xO [44 ] xl [36] - \ xO [60] x! [36, ' + 

We set T = 1 to avoid precision problems. 



T : = 1 

PulseSampling := T/8 



Taoie[ X ![i], {i, o, LaurontLK[L] [1] T / PulseSarapling -!}]]. 11 
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$ Aborted 



hessian = identityMatr ix [Length [var] ] ; 



Pulse[8] [0] [y_] 
Fulse[8] [1] [y_] 



x0[ Round [y / Pulse Sampling] ] ; 
xl [Round [y /PulseSampling] ] ; 



Temp Fun [x_] [y_] 



((AbsoluteValue[Pulse[8]] [x] [y] -1) A 2) //Expand; 



TerapFun2 [x_] [y_] 



= (AbsoluteValuelnter f earingregions [Pulse [8 ] ] [x] [y] ) ; 



TempFun[{-l, -1, -1, -1, -1. -D)[0.5T]; 



ConvertToBitSeq[Depth_] [Num„] : = 
Module [{bit, n) , 




We calculate the 



CostGenerator [TempFun„] : = 
Module [ {Cost} , 



Cost = 0; 
PulseValue = 0; 




mprun[#3 [PulseValue] U, PossSe.q] // Apply [Plus, #] & / / i, / 256 &; 



Table [CostUpdate, {i, 1, T / PulseSampling // Round} ] ; Cost]; 
AmplitudeCost = CostGenerator [TempFun] * PulseSampling; 
BERCost = CostGenerator [TempFun2] * PulseSampling; 

Save["AmplitudeCost.m", {BERCost, AmplitudeCost, PulseSampling, T 

In this section we find the bandwidth of the signal 
Derivative of the Sync function is represented by 



Clear [T, PulseSampling] 
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x J 
Cos[x] _ Sin[x] 

DSinc [0.0] : = 0 

DSinc[*_] /, x != o := COS W _ S"I*J 

<< AmplitudeCost .m; 
The bandwidth is given by 



Temp In to [m_ / n_] : = 

( Pulseslpli ns ) [°Sinc [ pulses ^ pling (x - m PulseSa^plin^) ] 

DSlnCf Pulsesa.npling (X " ° ^^^ampling) ] , { x, 0, 9 T) ] 
Templntl n_] : = 

( Pulsesa,pl ing ) N ^^«te [DSinc [ (x - » Puis a Sapling) ] 

DSinC [ Pulseslpling ( * " a ^^ampling) ] * 
{x, 0, 7T}] 

Table [Templnt0[0, n] , {n, 0, 5}] 

General: :stop : 

Further output of NIntegrate: :ncvb will be suppressed during this calculation. 
{1.27B79xlO\ 5. 06325, '10.6271, 9.78983, 4.51144, -2.28746} 
Pulse Sampling 

8 

ApproxBandWidth = Sum T (x0 fi + 11 - »nr^n ao y B i 

Bandwidth = sum[x0[i] x0[j] Templnt0[i, j], {i o 71^ M n 7 m 
Sum[xl t i] xl[J] Templntlfi, j] , {i/o,' 53> / « . 5, M} 3 ? ' ' >] 

Save[»BandWidthCost.m", {Bandwidth, T, } ] 
BandWidthCost = (ApproxBandWidth - 0.6) * 2; 
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BERWeight =0.3 

: £:!'- »— w«t«*i*t - B^e^t,. 

0.3 

TotalCost = BandWidthWeightBandWidthCost ♦ *mplitude W eight AmplitudeCost + 
BERWeight * BERCost; 

dir = Directory [] ; 

SetDirectoryf. ./new2satellite ] ; 
<< DavidonFletcherPowell -m 
SetDirectory[dir] ; 

L := 8 

DavidonFletcherPo W ell[TotalCo S t, var, hessian, initsol, 10], 
(0.0207507. (DavidonFletcherPov,ell'Private-alpha$39228-* 2,77153}} 
(0 . 0164912. (DavidonFletcherPowell • Private-alph a $39228-> 2 . 45972} } 
{0 . 0142842. (DavidonFletcherPowell • Private . a l P ha$39228- 1 . 88521} } 
{0 . 0126876. (DavidonFletcherPowell • Private ••lph«$3922B-. 2 . 08508) } 
{0 . 0117343. (DavidonFletcherPowell • Private >alpha$39228. 1 . 92086) ) 
(0 . 0112982. (DavidonFletcherPowell • Private ..lph.$39228-. 1 . 55191} ) 
{0 . 0110916, (DavidonFletcherPowell • Private . a lpha$39228-* 1 . 44356} } 
(0.010974. (DavidonFletch e rPowell>Private'alpha$39228- 1.63222}} 
(0.010886, { DavidonFl e tcherPow el l«Priv a te-a lp ha$39228- 2.1399}} 
(0 . 0108318, (DavidonFletcherPowell • Private>alpha$39228- 2 . 14114} } 
tempo [t_] = LaurentC[8] [0] [t] ; 



texnpl[t_] = LaurentC[8] [1] [t] 

initsoic 
join[Ta - - 

Table [tempi [i PulseSampl 

ListPlot[%278[[2]], PlotJoined -> Ime, PlotRange -> All) 
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0.1 



0.05 
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- Graphics - 



(^plituaecost, BERCost, BandWidthCost, App.oxBandWidth) /. Th read[var -> initsol] 

{0.0687008, 0.00375044, 0.0126982, 0.712686} 

(WltudeCost, BERCost, BandWidthCost, ApproxBandWidt*} /. Thread[var %278[[2]3] 
{0.0208114, 0.0Q872422, 0.00492789, 0.670199} 
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(ApproxBandWidth - 0.6) *2 /. Thread [var -> initsol] 
0.0126982 

(ApproxBandWidth - 0.6) *2 /. Thread [var -> *278[[2]]] 

0 . 00492789 

Thread[var -> %278[[2]]]] 

InterpolatingFunction[{{0, 8.875}}, <>} 

OptPulse[8]-[l] = IntBrpolationl Laurent LK [L] [1] T / Pulse Sampling - 1}] /. 

Table [{PulseSamplingi, xl[i]}, u, na^ie* l JL 

Thread[var -> %278[[2]]]] 

InterpolatingFunction[ { {0, 6.875}}, <>] 

tom3 . Kodulatcm^t^ ^oro^es - 2. Modularise . PptP-1..] , 

mterpolatingFunction^val : Input vaXue {^} lies outside the range of data 
in the interpolating function. Extrapolation will be used. 

i ,„™,f value 1^-1 lies outside the range of data 
InterpolatingFunction: -.dmval : Input value \ 32 J 

in the interpolating function. Extrapolation will be used. 

t T „n„t value (iiM lies outside the range of data 
InterpolatingFunction: : dmval : Input value { 1( . ] 

in the interpolating function. Extrapolation will be used. 
G |urther : ou? P ut of InterpolatingFunction: : d,nval will be suppressed during this calculation. 
I.i.tPlot[ { R.[t»3]. X m[ to m3]} //Transpose, P^oined -> True. AspectKatio -> 1] 




- Graphics - 

{Am pli t udeCost, «HCo.t. BandWidthCost, ApproxBandWidth, /. Tbreadtvar -> i-ltiol] 
{0.0687008, 0.00375044. 0.01.26982. 0.712686} 

{Anp li t udeCos t . BERCost, Ba»dWid th Cost, ApproxBandWidtH) /- Tnreadtvar -> «7.[[2n] 
{0.0208114. 0.00872422, 0.00492789, 0.670199} 
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Plot[M98[t], {t, 0 /7 T} # PlotRange -> All J 

InterpolatingFunction: :dmval : Inout valnp /g 07^,, . 

data in the interpolate funSSSn^^^ ° f 

InterpolatingFunction: :dmval : Inout valn^ /s onmoi 

data in the interpolate f unc^ion^ripo^t^n wif? ^ ""^ ° f 

InterpolatingFunction: rdmval : Input value /6 97*341 t *. • ^ 

data in the interpolating function^Extrip^lition win bfu'sfl. ^ ""^ ° f 

General: :stop *: 

Further output of Inter P olatin gFu nction : : a^val will De SU p preS sea d urin g this calculation 
0.2 



0.15 




- Graphics 



? Modulator 



? Interpolation 

^PP-x^ object which represents an 

{{xl, fl), {x 2, f2>? > or Ttt ff f^u The . data can have the forms 

xi are taken to have values 1, 2,'.T where ln the second . case, the 

Table [ 

Thread[var -> ?6247 [[2]]] 

Clear [Pulse] 

Sine [x_] /; x == 0 := 1 

Sine [x_] /; x != 0 := Sin[x]/(x) 

Pulse Sampling := t/8 

lullllt] M°} S " x °r R ound[y / Pulse Sampling] ]; 

Pulse[8] [1] [y_] := xl [Round [y/PulseSampling]]; 

ModulationValue[FiltPulse[81 1 rf-1 -1 t -i 

L J J Ll 1# 1 ' - 1 ' -1* -If "1/ -1}][0,5T] // N 

-0.754091 1 



ConvertToBitSeq[Depth__] [Num ] : = 
Module [{bit, n), 
n = Num; 

ConvertMext : = ( bit 
( ~ * " - 



*Z7 r** := (b±t = Mod t n ' 2]; n = Floor [n/ 21- bit> • 
Table [ConvertNext, (i, X, Dep th } ] // Reverse) 7^(0 3 ; 
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ConvertToBitSeq[8] [2] 
{-1, -1, -1. -1- - 1 ' - 1 ' 11 " 1} 

Table[i, <i, 1, 2*8}] // Map [Convert ToBitSeq[8] , #]&; 

Select [Table[i, (i. 1. 2-8}] // Map[ConvertToBit S eq[8] , #]*, (#[[5]] == 1)6] 

:i; :i: -i: 1: -i: d? i* i' i' "i 1 -!? 
:$; :i; :i: i: i: ^]'^-^^ -\\'- 

{-i, 
{-i. 

{-1, -3 
{-1. 
(-1. 
{-1, 
{-1. 
{-1. 
{- 
{- 
{- 
{■ 

{-: 
{- . 

iH : . l: l: I r. W^F**'" " " 

{1< -1, -1, -1. 1. 1. -1. -1). {1. - 1 ' "V j 1 '. 1 !* -1} {1. -1. -1. 1. 1. "I- "I- 

{i, -l, -l, -l. i. i. i. i>- t 1 .; - 1 : _1 ' x ' x : x . i it {1, -1. -1. 1. 1. -i}. 

1. -1, -1, 1. 1. -1. 1. I 1 *'* 1 ' 7 1 ' TVl 1 1 -1, -1,1. 1, 1. 1. 1}, 

i, -l, i, -i- 1. -i. -i. < x - -\' VV V A' A'} . ii -i. i. -i. 1. 1. - 1 - *>< 

\i, -i, l, -i. l. i . i. : i. 1. -1; i; 1; ^'i^Vi. -i. i. 1. i. -i t r 1 *' 

{1. 

[I] -i, i, i, i. i. .-i. -!>• tV "V 'i 1 '/' x' -i' -a.) i'. -i. - 1 ' - 1 ' 1} ' 

1, -1, 1, 1. 1. 1, 1. 1). {X. X, -X. ^ "i'.i'i 1 , {i. i, -1, -1, 1, 1. -1. -X}. 

{1, 1. 1. -1,1. X, 1, 1), {X, X 1. 1, 1,- . -J. ^ '{^1,1,1, 1,1, -1, -1}, 

• \l: 1:1: i: i: i^'-l: ^ < J - i - ^ ^ x - 1}} 
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xl = Map[ModulationValue[FiltPulse[8]] r#l TO 5 m & 

{ S:JJ"SI: o: 8 7 2 4 9 9 75' S'JS"2I' o-JJJ?il- S- 7 ^ 482 - °- 7 ^^ o. 7 n482. 0.711435 

0.776885 0 776839 6 776791 0 nnlila V2? 012 ' °- 792 °59< 0.792012 0.791965 
0-802178 O. 80213i; 8." 2o!oU' 0*ISf?I 8 ' 0 nlllY 1 ' 0 • 7 391- 0 . 739oi 4 ' 

0.759521. 0.759475 0.759428! o!759474 6 72^738 V^? V " T^ 393 ' °- 76 4347, 
0.819541, 0.819495, 0.819447 0 819494 07fl^^' ^'IIU*!' °- 72 1737, 0.72169, 
0.78714, 0.787094, 0 . 78704^ 0 787093 0 7««« 'n°^ If iS? 4 '«° " 781 - 756 ' 0-78171 
0.791922, 0.791876, 0.791829, 6 79^875 6 llAtk V«S?b£ VS? 355 ' °-749309 
0.749266, 0.74922, 0.749173, 0.74921 O 71 « 'n i? 5, « 0 ; 754138 ' 0.754091, 
0.829796, 0 82975 0 87<37nt n aoa^c,' „Z 11482 ' 0.711529, 0.711482, 0.711435 
0.776885 0 ??6M9, 6 "2?9i 6 7^6838 i^Pini 0 - 792 °l 2 ' 0.7919il,' 

0.802178, 0.802131 0.802084 0 2 'o 7 L 01 'A 0.7391, 0.739054, 

0.759521, 0.759475 0 75942R n " t^o^i n ' °- 764 44, 0.764393, 0.764347 

0.819541,' 0:819495! 0 bUH*' 0 Bllilt' o'vl^ll' Z^IU^' 0 ■ 0 .72169 , 

0.78714, 0.787094, 0.787047 0 787093 0 7««f 7 ' ft °4 7 fi28 4 ' °- 7817 56. 0.78171 
C791922, 0.791876. 0. 7,1*2* KXllh?^ 

Limit [D [Sin [x] /x, x], x -> 0] 



? Limit 

Li mit[ exp r . x - >x0] finds the liniting value of ^ ^ x approaches ^ 

««p[(l/2 - c/2Erf[±])*, xl] // A ppl y[Plus , # ]W/#/128t; 

BGr[a_] : = Evaluate 8] 
Ber[0.01] 

0.495 

Erf [infinity] 
1 

D[Erf [x], x] 



? Erf 

Erf{z] gives the error function erf(z) E rff 7 n 

function erf (zl) - erf(zO). erxt2 >- ErffzO, zl] gives the generalized error 
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Plot [Erf [x], {x, 0, 10}] 



2 




4 



6 



8 



10 



0.9999 I 
0.9998 
0.9997 
0.9996 

- Graphics - 

? FiltPulse 
Global * FiltPulse 
T := 3 / 812500 
Null 
?? T 

This is the symbol period 
?? T 

This is the symbol period 
T : = 3 / 812500 

Clear [T, x0, xl, x2, x3 # x4 , x5, x6, x7 ] 



* LaurentFunctions.nb 

) ; — 

BeginPackage [ "LaurentFunctions » « ] 
T::usage = "This is the symbol period" 
BT::usage = "This is the usual product" 
h::usag e = -This is the raw gaussian pulse" 
»: .-usage = "This denotes modulation index Pi" 
*::usage ="*[L,t] is Laurents function" 

hFiltered:: usage = "This is the filtered gaussian pulse" 

"^KSil-XS V^^^,^; ~^g«^. # -- -lowing 

S ; :usage = "S = Sin[3] » 
•7: :usage s " J - e i * w 
C: : usage = "C = Cos [$] " 
M: : usage = "13 = 2 L " 1 ■• 

fJodulationlndex: rusage = "Modulationlr.dex = h» 

Laurents: : usage = "Laurents [L] [a]'[t] = Sin[ *[L,t ♦ aTJ]/.s- 
r-aureatC: tusage = -Laureate [I.] [K] [t] is Laurents C K . t " 
AlphaKI::usage= "AlphaKI [LL] [K,i] is Laurents a*.* - 
LaurentLK: :usage = »LaurentLK[L] [K ] g ives the 6upport of Ck ^ . 
The start of Modulator Definitions 

ANXInitialStateSetUp: :usage = 

"^J^SSI^?^ -t. up the sequence 

constelation point sonified bv ™° ^ Went tn * ou 9 h »«t to the 

point specified by AccumulatedPhase which is really a„ 0 « 

A K N ::usage = " A K N[ L ]t K H{ State, A ccu 1 nulatedPhase )] defines A„ B in terms of A 0 .„ " 
ModulatingPulse: r usage = 

"The Pulse is assumed to have the following structure Pulse [L] [KJ [t] - 
NumberOfCurves: :usage = -The number of pulses used by the modulator" 
Samplinglnterval: : usage = 

"SamplingXnterval is the interval between samples of the output of the modulator- 

^^ir^^^^^^; ^ }m '^ of bits that are assumed 

StartingQuadrant:. - usage = "StartingQuadrant = A,._ ia nd is a number" 

Modulator: susage = "Modulator r LI rB-'t-c A ~ ^ 

ii^'^srrs :ss T u»\ 2 d^ h erofcu r es - ^^^^ 

he Pulse is assumed to have the following structure Pulse [L] [K ] [t] " 

Start of the Receiver Functions 
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FiltPul B e::usage="The default pulse and is called by FiltPulse [L] [K] [t] " 

Receiver: :usage = " Receiver [L] [InputSeq,OptsJ assumeB the following 

FiltPul S e,Syn=Sainple -> 0 , NunberOf Curves -> 2 

- -ppoeivorProperfl.] [ startingQuadrant, InitialState, 

Begin [ " * Private * w ] 
VLog[2] 



a : - 



2 7TET 



h[t_] := — 

V 2 7T a T 



Ideally we would and did define the effect of the convolution of h [t] by the formula below. However, this version of 

Mathematica gives an error. 

ri h[t - r] - . 

fcfii[t_] := Release [Module [{r>, J ^ <^JJ 

hFiltered[PulseWidth__] [t_] := Module [ {xl, x2, x3>, xl = 

rvhrtl-r] r PulseWidth PulseWidth T 
<N[#1, 40]*) [*«*!•[{«, J 2 2 20^' 

x2 = Interpolation [xl] ; x2 [t] ] 
S& :=N[ModulationIndex7r] 

C : = Cos [B] ; 
S := Sin [ft]; 
j := ( * // Chop) ; 
M : = 2 L ~ 1 ; 

PhaseAngle[L_J [t_J /; t S 0 : = 0 
PhaseAngle[L_] [t_] /; t s LT : = ff 

PhaseAngle[L_] [t_] := LT 
PhaseAngle[ LH t_] = «odule[{xl, x2, x3, x4, «5. x6), xl = hFiltered [3 X. «] [ tl - — ]; 

x 2= Table[{t2, . f Evaluate^!] dtl}, {«. 0, IT. ] ; Xnterpolation [x 2] [t] ] 

t_]/;0<t<LT:= PhaseAxigle[L] [t] 
t_] /; 2LT > t * LT s= « - PhaseAngle [L] [t - LI] 

We need to put this to avoid the function being extrapolated 

*[L.,t_]/;!(0<t<LT)«K2I.T>tiLT) :=0 
I,a»tentS[L_] [n_] [t_] : = Sin[ * [L, t + n T] ] / S 
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AlphaKI [LL_] [K_, i_] /. (0 < i < I>L) && (0 <= K < 2"- 1 ) : = 

Module [{xl, x2, x3, KNum} , xl := (x2 = Mod [KNiun, 2]; « ~ * 2 , x2 \. 

KNum = K; x3 = Table [xl, {ii, o, LI» - 1} ] ; x3 jfij J * 

LaurentLK [L_] [k_] : = 
Module[{xl}, *i = Table[ M 2 - AlphaKI [*] [K , ii] ) - ii, x, L - 1} ] • Min [xl] ] 

LaurentC[L_] [K_] [t_] /; 0 <= K < 2* : = 

LaurentS[L) [0] [t] f] LaurentS[L] [ii + L AlphaKI [L] [K, ii)] [t] 
ii.i - 

The start of the modulator function 

ANKInitialStatcSetUpCl,,] [K_] [XnitBitSe.L, AccumulatedPhase ] : = 

acuphase = AccunulatedPhase; 
UpdateSeq : = 

"acuphaseV^c^ (Vl. 1-1)], 

Table [ UpdateS eq . P {i,\, iS^SSS^^ 1 ]' iaitbitS<S « = Rest t ^itbitse,] ; xl'] ; 

AKN[L_) [K_] f{State_, AccumulatedPhase )] -= 
AccumulatedPhase - Sum [state [ [i ♦ 1] ] AlphaXl [X> ] [K, i] . {i , x, L - 1} ] 

Modulator [BitSeq_, Opts 1 : = 
Module [ J 

x^fsSiS^^alf/^ ™«. Pulse}, 

X3 a - Startf 1 o 1S ^ te '/ { ° PtS) OPtio^3 [Modulator]; 
" ? tar ^ n 30 uadr ant /. {Opts} /. Options [Modulator]; 

PulL m -T n - arVal/ - {0ptS) Options Modulator 
CurvL = M ^ ul « lB J»«l"/. {OPts} /. Options [Modulator], 

Acc^ulItid^asff^T 8 '* [«"**.*oi]' > - 1; 

seq = BitSeg; 
Table [ 

f ™^\;.^*^ «. 0, curves,]; 

AccumulatedPhase = AccumulatedPhase + First [seal- U ' 
seg = Rest[seg); 1 9J ' 

Table[ ^op?AK^^^^ 

x6[r_] = Sum [Sum [ ( a) *™«.t., K] tfl . a, , ^ + 

Table fxit;! 0 ' r LaU n en B f IjK[I ' ] [K] " 1}] ' < K ' °' Curves}] ; 
Table [x6 [r] , {t, 0, T - x4, x4}]] ; 

Table[x5, {kk, 1, Length [BitSeq] } ] //Flatten] 

ModuiaifngP^rsi 1 : FiltPuise'. ^S^l'/^ 1 o^ M ^i' r o^'c^; X :> 1 2 ) 1 , , ' SamPlinSlnte - 

Receiver [L_] [inputSe<a_, Opts ] : = 

Module [ {xl, x2, x3, x4, x5, x6}, 

*l = StartingQuadrant/. {Opts} /. Options [Receiver].; 
x2 = XmtxalState / . {Opts} / . Options [Receiver] ; 
X3 = Samplmginterval/. (Opts} /. Options [Receiver] • 
*4 = ModulatingPulse/. {Opts} /. Options [Receiver^/ 
«5 = SyncSample / . {Opts} /. Options [Receiver] ; 
x6 = NumberOf Curves /. {Opts} /. Options [Receiver] ; 
ReceiverProper[L] [xl, x2, x3, x4, x5, x6, InputSeq] ] 
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R.c.iverProper[X. ] [ StartingQuadrant., InitialState,, S^plinglnterval., 

^ividsTf'par"^ «D // Transpose ,/ # [ [SyncSa^le ♦ !])*; 

E ^dSlatIn e Pulse [L ] [0] [ (Laurent!* [I.] [0] / 2) T + SyncSa.ple Saznplinglnterva!] , 
.fl" = StartingQuadrant; 



Bcjn = 0; 
TD = {}; 

eeq = ReceivedSeq; 
ReceiveNext : = 
Module [{xl, x2}/ 



xl = ( <-l) son J* First [seq]) // Im; 
seq = Rest [seq] ; sgn = Modfsgn + 1, 2]]; 



seq = Rest [seq]; sgn - nou^y^ # ' i , ,, 
Table [ReceiveNext, {i, 1, Length [ReceivedSeq] }] ; 

ID] 
End[] 



EiidPac3tage[] 



